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ABSTRACT 

We have completed two years of photometric and spectroscopic monitoring of a large 
number of active galactic nuclei (AGNs) with very high accretion rates. In this paper, 
we report on the result of the second phase of the campaign, during 2013-2014, and the 
measurements of five new H/3 time lags out of eight monitored AGNs. All five objects 
were identified as super-Eddington accreting massive black holes (SEAMBHs). The 
highest measured accretion rates for the objects in this campaign are > 200, where 
^ = M./LEjdC -2 , M, is the mass accretion rates, L^d is the Eddington luminosity 
and c is the speed of light. We find that the H/3 time lags in SEAMBHs are significantly 
shorter than those measured in sub-Eddington AGNs, and the deviations increase with 
increasing accretion rates. Thus, the relationship between broad-line region size (R Hl3 ) 
and optical luminosity at 5100A, i? H/3 — L 5100 , requires accretion rate as an additional 
parameter. We propose that much of the effect may be due to the strong anisotropy 
of the emitted slim-disk radiation. Scaling R H6 by the gravitational radius of the black 
hole, we define a new radius-mass parameter ( Y ) and show that it saturates at a critical 
accretion rate of = 6 ~ 30, indicating a transition from thin to slim accretion disk 
and a saturated luminosity of the slim disks. The parameter Y is a very useful probe 
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for understanding the various types of accretion onto massive black holes. We briefly 
comment on implications to the general population of super-Eddington AGNs in the 
universe and applications to cosmology. 

Subject headings: accretion, accretion disks - galaxies: active - quasars: supermassive 
black holes 


1. Introduction 

Reverberation mapping (RM) experiments, which measure the delayed response of the broad 
emission line gas to the ionizing continuum in active galactic nuclei (AGNs), are the best way to 
map the gas distribution and derive several fundamental properties such as its emissivity-weighted 
radius and the black hole (BH) mass. The method was suggested by Bahcall et al. (1972) and 
its theoretical foundation explained in detailed in Blandford & McKee (1982). Numerous RM 
experiments, since the late 1980s (e.g., Clavel et al. 1991; Peterson et al. 1991, 1993; Maoz et al. 
1991; Wanders et al. 1993; Dietrich et al. 1993; Kaspi et al. 2000; Denney et al. 2010; Bentz et al. 
2009; Grier et al. 2012; Du et al. 2014; Wang et al. 2014a; Barth et al. 2015), have succeeded in 
mapping the emissivity distribution of H/3 and several other lines in the broad-line region (BLR) in 
more than 40 AGNs. BH virial mass measurements based on the method was shown to be consistent 
with the M, — cr* relationship, or with measurements based on stellar dynamics in local objects, 
where cr* is the stellar velocity dispersion in the host galaxy (Ferrarese et al. 2001; Onken et al. 
2004; see Kormendy & Ho 2013 for an extensive review). This method is an efficient and routine 
way to estimate BH mass at basically all redshift, at distances that are well beyond the resolving 
power of all ground-based telescopes. The RM experiments have been extensively discussed and 
reviewed in the literature (e.g., Peterson 1993; Kaspi et al. 2000 for earlier works and Shen et al. 
2015 and De Rosa et al. 2015 for more recent results). 

Most RM experiments, so far, have focused on the time lag of the broad H f3 line (r Hj3 in the 
rest frame) relative to the AGN continuum luminosity (XL\) at rest-frame wavelength of 5100A 
(hereafter Lsioo)- Results for 41 such measurements, based on AGN luminosity that is corrected for 
host galaxy contamination, are summarized in Bentz et al. (2013). They lead to a simple, highly 
significant correlation of the form, 

R Hf3 ~ a ltd , (1) 

where R H/3 = cr H/3 is the emissivity-weighted radius of the BLR and I 44 = T5ioo/10 44 erg s _1 . We 
refer to this type of relationship as the R Hj3 — L 5100 relationship. The constants a and /3 differ slightly 
from one study to the next, depending on the number of sources and their exact luminosity range. 
For the Bentz et al. (2013) work, a = 33.65 and /3 = 0.533. Scaling relationships based on other 
emission lines have been used to estimate BH mass in high-redshift AGNs. These are either based 
on the Mg nA2798 line, which is scaled to the H/3 line (e.g., McLure & Dunlop 2004; Vestergaard 
&: Peterson 2006; Shen et al. 2011; Trakhtenbrot &; Netzer 2012) or the C IVA1549 line, for which 
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few direct lag measurements are available (e.g., Kaspi et al. 2007). So far, there is not enough 
information about the dependence of the R Hft — £5100 relationship on accretion rate or Eddington 
ratio EBoi/^Eddi w h ere -^Boi is the bolometric luminosity, LecM = 1-5 x 1O 38 (M./M 0 )erg s _1 is the 
Eddington luminosity for a solar composition gas, and Af. the BH mass. 

We are conducting a large RM monitoring campaign targeting high-accretion rate AGNs. The 
aims are to understand better the physical mechanism powering these sources, the dependence 
of the R h/3 — T 5100 relationship on accretion rate, and the possibility to use such objects to infer 
cosmological distances especially at high z (Wang et al. 2013). We coin these sources “super- 
Eddington accreting massive black holes” (SEAMBHs). Earlier attempts in this direction, based 
on monitoring narrow-line Seyfert 1 galaxies (NLSls), failed mostly because of the small variability 
amplitude of the selected sources (Giannuzzo & Stirpe 1996; Giannuzzo et al. 1998; Shemmer & 
Netzer 2000; Klimek et al. 2004). Bentz et al. (2010) summarized results for a few NLSls, and 
found relatively small LBoi/^Edd for this class of objects and time lags that are not significantly 
different from the ones given by Equation (1). Our monitoring campaign started in 2012 and its 
first two phases are already completed. The ~ 20 targets observed so far are listed in Table 1. 
Results from the first phase have been published in Du et al. (2014; hereafter Paper I) and Wang 
et al. (2014a; hereafter Paper II). We also studied the time-dependent variability of the strong Fe 
11 lines in nine of the sources, and the results are reported in Hu et al. (2015, hereafter Paper III). 
We characterize SEAMBHs by their dimensionless accretion rate, = M % jLy^cr 1 , where Af. is 
the accretion rate. Typical values of this parameter range from a few to ~ 100 for the objects in 
the first phase of the project. Such high accretion rates are characteristics of slim accretion disks 
(Abramowicz et al. 1988) that are thought to power these objects (Szuszkiewicz et al. 1996; Wang 
Sz Zhou 1999; Wang et al. 1999; Mineshige et al. 2000). Paper II shows that such objects may 
eventually become new standard candles for cosmology. 

This paper reports the results of the second year of our RM campaign. Target selection, 
observation details and data reduction are described in §2. H/3 lags, BH mass and accretion rates 
are provided in §3. The new H/3 lags are discussed in §4, showing that for a given luminosity, the H/3 
lag gets shorter with increasing ^. In §5 we briefly discuss the implications of the new findings to 
the understanding of BLR physics and geometry, and to accretion physics. §6 gives a brief summary 
of the paper. Throughout this work we assume a standard cosmology with Hq = 67 km s _1 Mpc -1 , 
= 0.68 and Dm = 0.32 (Ade et al. 2014). 

2. Observations and data reduction 
2.1. Target Selection 

Unlike the first phase described in Papers I and II, in which most objects are NLSls, for the 
second phase we selected targets from the quasar sample of the Sloan Digital Sky Survey Data 
Release 7 (SDSS DR7) using the pipeline employed in Hu et al. (2008a). The objects have similar 
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spectroscopic characteristics to the NLSls from Papers I and II, in particular: 1) strong optical 
Fe ii lines, 2) narrow (< 2000km s -1 ) H/3 lines, 3) weak [O hi] lines (Osterbrock &; Pogge 1987; 
Boroson & Green 1992). The objects selected for the first year study are NLSls with extremely 
steep 2 — 10 keV continuum. For the second year sample reported here, we do not have X-ray 
data and use, instead, the dimensionless accretion rate , which can be estimated through the 
physics of thin accretion disks as formulated by Shakura & Sunyaev (1973, hereafter SS73). In 
such systems, the accretion rate can be directly calculated from the part of the spectrum where 
L u oc i^ 1 / 3 , regardless of the value of the BH spin (e.g., Collin et al. 2002; Davis & Laor 2011 and 
references therein). The only significant uncertainty in this estimate is the disk inclination to the 
line of sight, i (see more details in Paper II). We take cosi = 0.75 in this series of paper^j]. The 
standard thin disk equations give (see Paper II), 


A = 20.1 



3/2 

mf 2 , 


( 2 ) 


where mu = M./10'Mq. For thin accretion disks, we have LBoi/^Edd = where LboI is the 
disk bolometric luminosity, and rj is the mass-to-radiation conversion efficiency which depends on 
the BH spin. 

When comparing thin to slim accretion disks it is important to note that in both systems, 
the observed 5100A emission comes from large disk radii and thus is less influenced by the radial 
motion of the accretion flow compared with the regions closer in that emit the shorter wavelength 
photons. At these large radii, Keplerian rotation dominates, radiation cooling locally balances the 
release of gravitational energy through viscosity, and all effects arising from radial advection and 
the black hole spin can be neglected. As direct integration of the disk SED is not practical in almost 
all cases, due to the lack of far-UV observations, measuring L^oi/L^d directly is not possible, and 
Equation (2) is the best way to estimate which is directly related to the SS73 accretion disk 
model. 


To estimate the BH mass, we followed the standard approach and assume that the BLR gas 
is moving in Keplerian orbits and the rest frame time lags (r Hj3 ) provide reliable estimates of i? Hj3 . 
This gives, 


M. 


^hAfWHM 

L G 


1.95 x 10 6 / blr K 3 2 ti o Mq, 


( 3 ) 


where G is the gravitational constant, V 3 = Vp\VHM/10 3 km s -1 is the full-width-half-maximum 
(FWHM) of the H/3 line profile in units of 10 3 km s ^ 1 and = r H/ 3 /10days. The factor / BLR is 
calibrated from the known M, — a relation (e.g., Onken et al. 2004; Ho & Kim 2014). This is still 
a matter of some debate and the quoted uncertainties are large. The study of Ho &; Kim (2014) 
shows that for AGNs in host galaxies with pseudo-bulges, / BLR is smaller than in AGNs hosted by 


1 cosi = 0.75 represents a mean disk inclination for a type 1 AGN with a torus covering factor of about 0.5, 

assuming the torus axis is co-aligned with the disk axis. 
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classical bulges or ellipticals. However, / BLR gets larger for AGNs with higher Eddington ratios. 
It is not very clear what is the end result of these two opposite trends. For most SEAMBHs in 
Papers I and II, the host galaxies were observed by the Hubble Space Telescope ( HST ) and show 
indications for pseudo-bulges. They also have very high Lb oi/LEdd- For the present sample we do 
not have HST images but these sources also are of high Lb 0 i/Lecm- Given all this, we use / BLR = 1, 
as in Papers I and II, but note the large uncertainty on this number. 

Employing Equations (1) and (3), we estimated ^ from Equation (2) for quasars in the 
SDSS DR7 and selected sources with the highest ^. We discarded radio-loud objects based on 
the reported FIRST observations. We also chose a redshift range of z = 0.1 — 0.3 and L 5100 < 
10^ 5erg s ' 1 to make sure that the lag, as estimated from Equation (1), can be measured in 5 — 6 
months of monitoring campaign and magnitude r' < 18.0 to ensure high enough signal-to-noise 
(S/N). Details of the sources are given in Table 1. 


2.2. Photometry and Spectroscopy 

The second year observations described here are essentially identical to those of the first year, 
and readers are referred to Papers I and II for detailed information about the telescope and spec¬ 
trograph. In short, the telescope is a 2.4 m alt-azimuth mounted Ritchey-Chretien located in 
Lijiang, Yunnan Province, and operated by the Yunnan Observatories. We use the Yunnan Faint 
Object Spectrograph and Camera (YFOSC) with a back-illuminated 2048x4608 pixel CCD cov¬ 
ering a field of 10' x 10'. The differences between the first and 2nd year observations are: 1) we 
adopted a 5 w -wide slit, compared with the previously used 2.5" slit, to minimize the influence of 
atmospheric differential refraction, and used Grisrn 3 with spectral resolution of 2.9A/pixel and 
wavelength coverage of 3800-9000A. 2) Instead of the Johnson V filter used in 2012, we used SDSS 
r'-band filter to avoid the potential contamination by emission lines such as H /3 and Ha. 3) The 
observations now include photometry from the Wise Observatory, lrn telescope in Israel, where 
we used a Princeton Instruments CCD camera, an SDSS r '—band filter and a typical exposure 
time of 900 sec. The reduction of the Wise and Yunnan Observatory photometry data was done 
in a standard way using IRAF routines. For the Lijiang data, the reduction, and the procedures 
to measure 5100A and F(H/3), are given in Papers I while for Wise data, the flux measurements 
were done using point-spread function photometry. The light curves were produced by comparing 
the instrumental magnitudes to those of constant stars in the field (see, e.g., Netzer et al. 1996, 
for details). The uncertainties on the photometric measurements include the fluctuations due to 
photon statistics and the scatter in the measurement of the stars used. Detailed information of 
observations, such as, comparison stars etc. is provided in Table 1. We also list the critical S/N 
used as a low threshold to reject poor weather condition observations. 

The 5100A light curves used in the analysis were obtained by inter-calibration of the 5100A 
spectroscopy, r '—band Lijiang 2.4-rn photometry, and r '—band Wise Observatory photometry. We 
used the Bayesian algorithm described in Li et al. (2014), who assume that the AGN optical contin- 
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uum variations follow a damped random walk model (Kelly et al. 2009; Li et al. 2013). We applied 
a multiplicative scale factor and an additive flux adjustment to bring the different measurements 
to a common flux scale. The algorithm performs inter-calibration on all the datasets simultane¬ 
ously. This enables us to relax the requirements for the sampling rates and retain the highest 
achievable temporal resolution. We selected the Lijiang spectroscopy as the common scale. The 
best parameters for the inter calibration are determined by a Markov chain Monte Carlo (MCMC) 
implementation. In that way we get the inter-calibrated continuum light curves of photometry 
(Lijiang and Wise data) and spectral continuum (Lijiang data). Finally, the combined light curves 
are obtained by averaging all the points of a same day in the inter-calibrated light curves. All 
the continuum and H/3 light curves for the eight objects are listed in Tables 2 — 5, and shown in 
Figure 1. We also calculated the mean and RMS (root mean square) spectra and present them in 
Appendix A where we briefly discuss these data. 


2.3. Host Galaxies 

Most of the objects described in Papers I and II had previous HST images that were used to 
subtract the host galaxy contribution at 5100A. Such information is not available for the sources 
reported in this paper. Instead, we used the empirical relation suggested by Shen et al. (2011) 
based on SDSS spectra: L^oo/^ioo 1 = 0-8052 — 1.5502x + 0.912x 2 — 0.1577x 3 , for x < 1.053, where 
x = log (L 5 °Q 0 / 10 44 erg s _1 ) and LI^qq is the total emissions of AGNs and their host at 5100A. For 
x > 1.053 we have Lj?°Qo -C q, and the host contamination can be neglected. These authors 
obtained the geometric mean composite spectra for objects binned in logZ /5100 normalized at 3000A 
and found spectral flattening toward long wavelengths. The flatten part of the spectra depends 
on host contaminations. This empirical relation was obtained by the difference of the flatten part 
through comparisons among the spectra for quasars in 5100A luminosity bins. 

The above procedure was developed for the 3" SDSS fiber noted as x 3 „. This is considerably 
smaller than the effective aperture (5") used by us. Denote x 5 „ as the flux measured from our 
mean spectra and Ax = \x 3 „ — x „\. We find that on average Ax/x 3 „ < 0.3. We also find that 
the differences in L^qq/L^^ between the 3" fiber and our 5" aperture is less than 10%. The 
fraction is estimated by comparing Fy and Fy fluxes at 5100A in our sample, where Fy and Fy 
are measured from Lijiang and SDSS observations, respectively. In general, more than 90% of the 
light is contained within the 3" fiber for z ~ 0.2 quasars. Thus, the present subtraction of host 
contamination is quite robust, and the added uncertainty due to this difference is not larger than 
the intrinsic uncertainty in Shen et al.’s expression. 
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3. H/3 Lags, Black Hole Mass and Accretion Rates 

We used a standard cross-correlation analysis to determine the lags of the H/3 line relative to the 
combined 5100A continuum. The procedure is quite standard and was described, in detail, in Papers 

1 and II, where we also list all the relevant references. The uncertainties on the lags are determined 
through the “flux randomization/random subset sampling” method (RS/RSS) (Peterson et al. 1998, 
2004), the cross-correlation centroid distribution (CCCD) and cross-correlation peak distribution 
(CCPD) generated by RS/RSS method (Maoz & Netzer 1989; Peterson et al. 1998, 2004; Denney 
et al. 2006, 2010 and reference therein), which are shown in Figure 1. For a successful detection of 
H/3 lag we require: 1) non-zero lag from the CCF peak and 2) a maximum correlation coefficient 
larger than 0.5. 

Five of the eight sources show well-determined lags. All these sources are listed in Table 6 . The 
other three failed our above criteria. Figure 1 shows that the CCFs of two of them (J094422 and 
J100055) show two comparable peaks and no significant H/3 lags. The third source, J080131, has an 
unusual combination of line and continuum light curves. The continuum light curve shows a well- 
determined, broad feature followed by two major short ’’bursts” centered at around JD24566000+, 
120 and 140. The H/3 light curve shows a clear response to the first continuum dip and burst 
feature, but no response to the two shorter duration features. The CCF of the entire campaign 
has a low peak, at ?’ max — 0.5, which is consistent with a zero lag. Carrying the analysis for the 
first 70 days only, as shown in the second diagram of J080131 panel in Figure 1, we found a very 
significant peak with a centroid lag of 11 . 5+30 days (with a very high coefficient of r max = 0.81) in 
rest-frame. 

Our search of the literature shows that the unusual combination of line and continuum light 
curves observed in J080131 is rare [see similar, but not identical, behavior in NGC 7469 (Peterson 
et al. 2014) and in the UV light curve of NGC 5548 (De Rosa et al. 2015)]. We can think of 
various scenarios that could cause such an event in a SEAMBH with extremely high accretion rate. 
For example, the inner part of slim disks can be so thick that self-shadowing effects lead to strong 
anisotropy of the emitted radiation. The signal received and measured by a remote observer can 
differ substantially from the ionizing radiation that reaches the H/3-emitting clouds (see discussions 
by Wang et al. 2014b). Given the rarity of such objects, we have no way to test this idea and thus 
do not include this object in the remaining analysis. 

We used Equations 2 and 3 to calculate accretion rates, ^ and BH masses for the five sources 
listed in Table 6 . The details of the procedure are given in papers I and II and the distributions 
of ./# in the two groups of the previously and currently monitored objects are shown in Figure 
2. As clearly shown in the diagram, our Lijiang-observed sources occupy the high tail of the ^ 
distribution. This is not surprising given the way we selected our targets. We also show in Figure 

2 the equivalent width (EW) distribution of H/3. On average, our high accretion rate sources have 
lower mean EW(H/3). An anti-correlation between EW(H/3) and Lboi/LecM has been found in other 
samples with less accurate rates, in which BH mass is estimated by the i ? H/3 — L 5100 relation (e.g., 



Netzer et al. 2004). 


In Paper II we classified SEAMBHs as those objects with m = r\.M > 0.1. This is based on 
the idea that beyond this value, the accretion disk becomes slim and the radiation efficiency is 
reduced due to photon trapping and other effects. Since we cannot observe the entire SED, we 
have no direct way to measure LBoi/TEddi an ^ this criterion is used as an approximate tool for 
identifying SEAMBH candidates. To be on the conservative side, we chose the lowest possible 
efficiency, rj = 0.038 (retrograde disk with a = —1, see Bardeen et al. 1972). Thus SEAMBHs are 
objects with = 2.63. For simplicity, in this paper we use = 3 as the required minimum. 

Later on, in §5, we introduce an empirical way, based on our own mass measurements, to define 
this group more accurately. 

Table 6 lists all observables, Af, and .M for the first (SEAMBH2012) and second (SEAMBH2013) 
year observations. All objects, except for MCG+06 —26 —012 (see note in the table caption), show 
> 3, indicating they are SEAMBHs. In Paper III, more accurate measurements of H/3 lags have 
been derived by using a scheme of simultaneously fitting spectra of host and AGNs. The newly 
measured lags are the ones listed table 6. Comparing with the previous measurements in Papers I 
and II, the updated lags are consistent with the previous ones, but the error bars are much smaller. 
Including our new observations, the majority of SEAMBHs with RM-based mass measurements 
come from our Lijiang RM campaign. 


4. H/3 Time Lags in SEAMBHs 

To understand better the fundamental differences between high and low ../M sources, we must 
use the measured correlation between time lag and source luminosity for all RM AGNs. For this, 
we have to take into account that some of the sources have been mapped more than once and hence 
there are various measured values for their luminosity, time lag, and FWHM(H/3), for the same BH 
mass. There are two possible ways to address this issue. The first is to average the BH mass from 
the individual campaigns, and then obtain mean (or median) values for L 5100 , L H(3 and the time 
lag. This scheme, which was used in Kaspi et al (2005), Bentz et al (2009) and other papers, is 
referred here as the “average scheme”. The second is to treat the independent RM campaigns of a 
single object as different objects (e.g., Bentz et al. 2013). We call this the “direct scheme”. 

The goal of the present paper is to find the basic properties of the central power house, 
presumably an accretion disk, as a function of BH mass and normalized accretion rate. The 
measured H/3 time lag is the tool we use to study these properties. This is changing with the 
ionizing continuum but it does not represent real changes in the gas distribution in the BLR during 
one observing season. Since the BH mass does not change on a short time scale, and the time scale 
to change the global accretion rate through the disk is also very long, we must be seeing short 
time scale fluctuations that do not change these fundamental properties but, nevertheless, affect 
the ionization of the BLR. To understand the population properties, we must give equal weights 
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to all objects, i.e., use the average scheme approach. This means that the mean luminosity is used 
to derive -M. This can be done during one campaign that lasts a few months, or during several 
campaigns that last a few years. The only exception is, perhaps, when the total time exceeds the 
dynamical time scale of the BLR. 

The direct scheme can be useful to derive other properties, such as the ionization distribution 
across a “typical BLR”. While this is not the goal of the present work, we have calculated all the 
fundamental correlations discussed below using this approach, too, and show them in Appendix B. 
We note that despite the fundamentally different approaches, which result in a different BH mass 
distribution (since several masses are counted more than once), the results are not very different. 


4.1. Continuum luminosity vs. BLR Size 

Table 7 lists the objects summarized by Bentz et al. (2013), plus a few newly mapped sources. 
Since some of the procedures used in earlier campaigns differ in several aspects, we decided to apply 
the same uniform procedure for all the sources. For objects with multiple RM measurements, our 
procedure for obtaining the average values is as follows: We first calculate M, for each of the 
individual campaigns and then average all these values, using weighted standard deviations, to 
get a mean M, and its error. This value is used to obtain individual values of ^ for the various 
campaigns, using L 5100 from this campaign. The mean for the object is the average of the 
individual and its errors include the standard deviation. For the mean luminosity and r H ^, we 
used logarithmic averages. 

Regarding BH mass calculations in individual campaigns, this is more problematic since some 
works used the RMS spectrum (e.g., Peterson et al. 2004; Bentz et al. 2009b; Denney et al. 
2010; Grier et al. 2012) while others prefer the use of FWHM(H/3) (e.g., Kaspi et al. 2005). Our 
procedure is based on the FWHM method (see explanation in Papers I and II). As shown recently 
in Woo et al. (2015), the scatter in the scaling parameter (/blr) derived in this method is very 
similar to the scatter in the method based on the RMS spectrum. 

We have therefore collected from the literature all the required information about the H (j line 
and carried out our own mass measurements for the sources where earlier mass measurements were 
based on the RMS method. In most of the cases, the agreement between the older RMS method 
and the value obtained here is within the uncertainty on the BH mass. In a few cases (five objects: 
NGC 4051, NGC 4151, PG 0804+761, PG 1613+658 and PG 2130+099) the deviation is larger 
which we interpret as unrealistically small uncertainty on the mass. The deviation is generally 
smaller than ~ 0.5 dex except for NGC 4051 where the deviation is ~ 0.7dex. Some more details 
for FWHM measurements of a few objects are given in the caption of Table 7. We list all the 
measured values and their averages in this tablco 


2 Following Bentz et al. (2013), we included in the uncertainties of the 5100A luminosities also the uncertainties 
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All correlations shown in this paper were calculated with the FITEXY method in the version 
adopted by Tremaine et al. (2002), where scatter is allowed for by increasing the uncertainties in 
small steps until y 2 is about unity (this is typical for many of our correlations). We also used the 
BCES method (Akristas & Bershady 1996) but prefer not to use its results since it is known to 
give unreliable results in samples containing a few outliers (there are a few objects with quite large 
uncertainties of .JK ). 

Figure 3a shows H/3 lags versus T 5100 for all the observations of mapped AGNs (i.e., multiple 
measurements for each object when available). Our objects are marked with blue (first year, data 
reported in Papers I and II) and red (second year, this paper). Using the data in Table 7 we get: 


log (^ H/3 / ltd ) 


' (1.45 ± 0.03) + (0.50 ± 0.03) log I44 

(entire sample), 

(4a) 

< (1.54 ±0.03)+ (0.53 ±0.03) log Z44 

(,A < 3), 

(4b) 

k (1.30 ± 0.05) + (0.54 ± 0.06) log I44 

{,J, > 3), 

(4c) 


with intrinsic scatters = (0.21,0.15,0.24) for (4a, 4b, 4c), respectively. The regression for the 
< 3 sample is almost identical to Bentz et al. (2013). Equations (4b) and (4c) show different 
intercepts for the Jt < 3 and the > 3 samples. Most of the SEAMBHs from our campaign are 
located below the i ? Hj3 — T5100 relation, increasing the scatter in the relationships. For more extreme 
SEAMBHs with Jt > 10, the regression shows log(i? H/ 3 /ltd) = (1.18 ± 0.07) + (0.48 ± 0.10) log ^ 44 , 
indicating they deviate from the < 3 group more than the Jt > 3 group. It is clear that the 
SEAMBH sources increase the scatter considerably, especially over the limited luminosity range 
occupied by the new sources. 

In the following analysis, we use the averaged lags and the averaged mass as listed in Table 
7 (note that unlike in paper-II, here we include in the analysis the few radio-loud AGNs with 
measured time lags and BH mass). In Figure 3c, we divided the population into two sub-groups, 
those with >3 (11 from earlier studies and 13 from our study) and those with Jt < 3. The 
diagram emphasizes that much of the intrinsic scatter is caused by a systematic deviation of the 
large ./# sources towards shorter lags. The regression calculations give 


lo g (^W^d) 


' (1.45 ± 0.04) + (0.50 ± 0.03) log I44 

(entire sample), 

(5a) 

< (1.55 ±0.04)+ (0.53 ±0.04) log Z 44 

{,A < 3), 

(5b) 

k (1.32 ± 0.06) + (0.53 ± 0.06) log U 4 

(,A > 3), 

(5c) 


with intrinsic scatters of = (0.21,0.16,0.22) for (5a, 5b, 5c), respectively. The slope of the 
correlation for the SEAMBH sample, those with Jt > 3, is comparable to that of sub-Eddington 
AGNs, but the normalisation is significantly different. Clearly, the H/3 region of SEAMBHs is 
smaller than sub-Eddington AGNs (see the different intercepts in Equations 5b and 5c). 


due to distances. 
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A more complete approach is to analyze the dependence of r Hj3 on L 5100 and together, 
however, several reasons prevent us from performing such an analysis. First, the current sample is 
still small; and second the luminosity range of SEAMBHs is very narrow. There are only 4 lumi¬ 
nous SEAMBHs (PG 0026+129, PG 1211+143, PG 1226+023, PG 1700+518) identified from the 
previous campaigns, but their = (4.5,6.9,5.0,12.0), respectively, are smaller than the average 
of this group. Thus, the search for a new relationship of the type of R H/3 = -R H/3 (T5ioo> will 
have to wait for future observations of luminous, high accretion rate sources. 


4.2. H (3 luminosity vs. BLR Size 

Figure 3c shows the new R Hj3 — L H/3 relationship. Such relationships have been used in the 
past to supplement the continuum-based relationships and to correlate the BLR size more closely 
with the ionizing continuum radiation (see e.g., Kaspi et al. 2005). The regression analysis gives, 

(1.30 ± 0.03) + (0.51 ± 0.03) log Lh +42 (entire sample), (6a) 

(1.36 + 0 . 04 ) + ( 0.52 ± 0 . 04 ) log Lh+42 {-4 < 3), (6b) 

(1.21 ±0.06)+ (0.53 ±0.06) log L h +42 (^>3), (6c) 

where ±h +42 = /10 42 erg s -1 , and <r ill = (0.20,0.18,0.22) for (6a,6b,6c), respectively. As in 

the continuum luminosity case, the inferred BLR size for SEAMBHs is smaller than the size of the 

sources with < 3 but the differences are smaller in this case. 

The R h/3 — relationship has been examined by Wu et al. (2004) and Kaspi et al. (2005). 
A comparison with the results of these papers shows similar intercepts but different slopes, with 
the present slopes being significantly flatter than 0.685 ± 0.106 found in Wu et al. (2004) and 
0.694 + 0.064 found by Kaspi et al. (2005). It is not at all clear that these are significant variations 
since the present sample is much larger and the distribution of sources along the L H/3 axis is quite 
different. The correlations involving L 5100 and L H/3 (Figure 3 and Figure 9 in Appendix B) show 
a very similar scatter but, for a given H f3 luminosity, SEAMBHs have shorter lags compared with 
low-accretion rate AGNs. Such differences may be related to the properties of slim accretion disks 
(Wang et al. 2014b), or other types of anisotropies (e.g., Netzer 1987; O’Brien et al. 1994; Goad 
& Wanders 1996; Ferland et al. 2009). 



4.3. 


'—dependent BLR Size 


To test the dependence of the BLR size on accretion rates, we define a new parameter, Ai? Hj3 = 
log (R H/3 /R H(3 r _ l ) that specifies the deviations of individual objects from the R H/3 — L 5 ioo/± H(3 
relationship of the sub-sample of < 3.0 sources (i.e., R H/3 R _ L as given by Equations 4b, 5b 

and 6b). The scatter of AR H/3 is calculated by a = |")>± (ARh+i — (ARh/?)) 2 /ivl , where N is 
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the number of objects and (Ai?n/?) is the averaged value. Figure 3 provides the value of a for 
comparison. 

Figure 4 shows plots of A R H0 versus . and Ai? H/3 distributions for the ,/0t > 3 and < 3 
samples. A Kolmogorov-Smirnov (KS) test comparing the two shows that the probability of the 
same parent distributions is p KS = 0.0054 for the R Hj3 — L 5100 case and p KS = 0.014 for R H0 — L H/3 
case. This provides a strong indication that the main cause of deviation from the old R H0 — T 5100 
relationship (or the R H0 — L 5100 relationship for the sub-Eddington AGNs) is the extreme accretion 
rate. Thus, a single i? Hj3 — L 5100 relationship for all AGNs is a poor approximation for a more 
complex situation where both the luminosity and the accretion rate determine this relationship. 
From the regression, we get the dependence of the deviations of R H0 from the i? H/3 — T 5100 relation 
in Figure 4a 

A R k0 = (0.36 ± 0.13) - (0.46 ± 0.10) log A (for A > 3), (7) 

with (Ji n = 0.07. From the R H/3 — relation in Figure 4c, we get 

A R h0 = (0.33 ± 0.15) - (0.38 ± 0.10) log A (for A > 3), (8) 

with ai n = 0.13. We have tested the above correlations also for < 3. The FITEXY regressions 
give slopes of around zero with very large uncertainties, A R H0 oc ^-°- 085 ± 0 051 an d A R H0 oc 
^-0.059±0.054 £ or p£g ure 4 a anc | 4 Cj respectively, implying that Al?, Hfl does not correlate with 
for ^ < 3 group. All this confirms that is an additional parameter controlling the R HPi — T 5100 
relation in AGNs with high accretion rates. 

To summarize, the new SEAMBHs observed in 2012 and 2013 significantly increase the scatter 
of the — L 5100 relation. We find that SEAMBHs have significantly shorter lags than those 
of sub-Eddington AGNs with similar 5100A luminosity, and the shortened lags increase with the 
dimensionless accretion rate. Given this, we recommend to use Equation (4b) or (5b) when trying 
to estimate the BLR size in sub-Eddington AGNs. We suggest that the dependence of i? H/3 on the 
dimensionless accretion rate could be a consequence of the anisotropic radiation of slim disks. We 
come back to this issue in the following section. 


5. Discussions 

The new results presented here illustrate that the small scatter in the earlier R H;3 — T 5100 
relationships (Kaspi et al. 2005; Bentz et al. 2013 and references therein), only applies to low 
accretion rate AGNs, which, according to our definition, are sources with ^ < 3. The addition 
of high accretion rate sources, those referred to as SEAMBHs, changes this picture. The scatter 
increases and there is a clear, statistically significant tendency for AGNs with higher accretion rates 
to show smaller i? H/3 . Moreover, the deviation increases with ^. This holds for both T 5100 and 
L H/3 although the deviations are larger when using the 5100A luminosity. The deviations are most 
noticeable over the luminosity range where most SEAMBHs are concentrated, between 10 43 and 
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few x 10 44 erg s _1 . Below we explore the physical implications of the new results focusing on the 
earlier suggestion (Paper II) that SEAMBHs are powered by slim accretion disks with properties 
that are considerably different from those of thin disks. 


5.1. and the Origin of the Shorter Time Lag 

The deviations from the earlier R Ha — L 5100 relationship raise several interesting possibilities 
that can explain the new results. The first possibility is that the ionizing flux reaching the BLR, 
which determines its RM-measured size, is different from the flux inferred by a remote observer due 
to an unusual SED and/or the anisotropy of the radiation emitted from the disk. This possibility 
has been mentioned in the past with regards to the covering factor and the level of ionization of the 
BLR. Netzer (1987) presented a detailed model where the combination of an angle-dependent SED 
of geometrically thin accretion disks (see also Czerny & Elvis 1987; Li et al. 2010 for geometrically 
thick parts of the disks), with an isotropic X-ray source, results in a weaker ionizing continuum 
along the plane of the disk compared with the pole-on direction. This will result in a complicated 
angle-dependent level of ionization. The recent work by Wang et al. (2014b) suggests that large 
SED variations, especially in the ionizing luminosity as a function of the polar angle, will be most 
noticeable in slim accretion disks where the anisotropy of the disk is more extreme because of 
self-shadowing effects. Such a geometry would lead to a smaller emissivity-weighted BLR radius 
at large angles relative to the polar axis. The second possibility is that much of the effect is due 
to the change of the bolometric luminosity with polar angle due to the reasons discussed above. 
This will reduce the dust sublimation radius that has a large effect on the measured RM size (see 
below). Obviously, the covering factor of the BLR, which determines L H/3 , and the anisotropy of 
the Balrner line emission, could also affect these relationships. 

The search for a direct ionizing luminosity indicator initiated several attempts to replace L 5100 
by L H0 in the R Hj3 — L 5100 relationship because H f3 responds more directly to the ionizing luminosity 
(e.g., Peterson 1993). Kaspi et al. (2005) show that this substitution had no affect on the deduced 
BLR size. Similar ideas have been proposed by Wang and Zhang (2003), Wu et al. (2004), and 
others. Our new observations hint to the possibility that the changes in the L H/3 relationship relative 
to the relationship for the ^ < 3 group are smaller than in the L 5100 case. Unfortunately, the 
relatively small size of the SEAMBH sample prevents us from quantifying these ideas. 


5.2. Accretion Disks and the Maximum Size of the BLR 

5.2.1. Geometrically thin disks 

Here we explore the possibility that the observed changes in BLR as a function of are due 
to the transition from thin (./# < 3) to slim (^# > 3) disks. To start, we express the measured H /3 
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lags in units of the gravitational radius, R g = 1.5 x 10 5 (Af # /M Q ) cm. Since R H/3 oc l b Ai (Equation 
1), and ^44 oc (Equation 2, recast), we expect 

r Hj3 = ^ = 1.9 x 10 4 .# a 35 ?nf a29 , (9) 

where we use Equation (5b). Equation (9) concisely mergers the empirical R H/3 — L 5W0 relation 
with the basic expectations of geometrically thin disks. 

We now defined a new parameter, Y = m 7 29 r H/3 , which we coin “the radius-mass parameter”, 
i.e. 

Y = 1.9 x 10V/°' 35 . (10) 

We calculated this combined parameter for all sources in order to test whether it depends on 
accretion rate. This test would provide clues about the various modes of accretion and their 
relationship with the size of the BLR. 

Similar to the definition of Ai ? Hj3 in Section 4.3, we define AX = log (X/X ieg ), where X is 
the measured value for individual sources and A reg is the value obtained for the group from the 
regression analysis. In this section we consider two cases, X = r H/3 and X = Y. The various tests 
performed on the accretion rate groups are shown in Figure 5 and discussed below. 

Our first test for r H/3 is performed on the two accretion rate groups, ^ < 3 and > 3. Figure 
5a shows the r H;3 — relation. For .Yt < 3 objects, we find r Hj3 oc ^°- 50 ± 0 - 08 with < 7 i n = 0.24, 
which is steeper than expected (Equation 9). By contrast, r Hi3 oc ^o. 35 ±o.i 2 for the Yt > 3 group, 
consistent with Equation (9). It is clear that the r Hj3 — Yt relation is markedly different for the 
two groups, suggesting that SEAMBHs are physically distinct from sub-Eddington AGNs. Ar H;3 
and its distributions are shown in the lower panels of Figure 5a. A KS test gives the probability of 
p KS = 1.3 x 1CT 3 and demonstrates the differences again for the > 3 and the < 3 groups. We 
suspect that the deviations from the expected dependences on ^# 0 35 are due to the dependence 
on 777,7 in our sample that contain a large range of BH mass. 

Figure 5 b shows the Y — relation, which has much smaller scatter than the r H/3 — ./# rela¬ 
tion, for both groups. The regression gives Y oc ^#°- 41±0 - 07 for the < 3 group with cr m = 0.15, 
which, within the errors, agrees well with Equation (10) confirming the idea of thin accretion disks 
in objects with ..-W, < 3. The diagram clearly shows a change of slope at a transition accretion rate 
of Jlc > 3 from a rising to an almost horizontal constant (saturated) value. This transition corre¬ 
sponds to a saturated— Y of y sa t ~ 3 x 10. This is consistent with the idea that beyond this critical 
rate, which we name ^# c , the central power-house is no longer a thin accretion disk. Obviously, 
such a test could not be performed without the new observational data from SEAMBH2012 and 
SEAMBH2013. This is not entirely unexpected since previous works specifies in the introduction 
to this paper, and explain in more detail below, show that the SS73 disk model becomes invalid 
above a certain accretion rate (e.g., Laor & Netzer 1989). 
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5.2.2. The Transition Accretion Rate 


Figure 6 shows the determination of the transition accretion rate and the saturated— Y. To 
determine from the data, we first test the dependence of the index b of the R H/3 — T 5100 relation 
on ^# c , and find b = (0.527 ± 0.038,0.528 ± 0.030,0.533 ± 0.028) for samples of objects with 
^ c = (3,10,20). These changes are so small that we decided to adopt b = 0.53 for all the 
groups. In order to obtain .M c and Y^ a t, we use the following functional fornjfl. 


Y = 



< 


;> .^c), 


( 11 ) 


where Igat, ko and are to be determined by the data. We make use of the Levenberg-Marquardt 
method to fit the entire sample with the inclusion of error on jjt through Monte-Carlo simulations 
(e.g., Press et al. 1992). In the fitting, we increase the uncertainties of Y until y 2 = 1. This gives 

Jt c = 13.8+g 9 ! 6 , y sat = 3.3ig;jj x 10 3 4 * , and k 0 = 0.27 ± 0.04. (12) 

The value of k 0 is consistent, within 2cr, with the expected value of 0.35, providing empirical 
support, from the basis of RM experiments, for the geometrically thin disk model. 

The analysis shows high uncertainty on the value of This is due to the limited numbers 
of objects with > 13.8 (13 objects). The precision we can determine the values of .Y£ c and y sat 
would most likely improve when large samples of high accretion rates become available. 


5.2.3. Slim Accretion Disks 

To understand the dependence of the radius-mass parameter on the accretion rates, we need 
to consider the properties of slim accretion disks. Such disks are the result of the increasing mass 
accretion flow through a thin accretion disk which leads to a build-up of radiation pressure close to 
the BH, the thickening of the disk, and the trapping of the photons that are trying to escape from 
these regions. Such systems have been discussed in numerous papers (e.g., Abramowicz et al. 1988; 
Wang Sz Zhou 1999; Mineshige et al. 2000; Watarai et al. 2001; Sadowski 2009, and references 
therein). The general result obtained from these studies is a simple dependence of the bolometric 
luminosity of the system linearly on the BH mass and only a weak logarithmic dependence on the 
mass accretion rate (so called “saturated luminosity”). The characteristic saturated bolometric 


3 We also tried to fit the data with a broken power-law, Y = Y sa ,t^ kl [l + (.-^'/^ c ) k2 ]~ 1 , where k 1 , k%, Rat and 

can be obtained from the fit. This did not result in a satisfactory fit. The reasons are: 1. Many more objects 

with low .-ft regardless of the exact value of ./M c . 2. The overall range in for low accretion rate sources (more than 

3.5 dex) is much larger than for high accretion rate sources (less than 1.5 dex). 
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luminosity can be expressed as (Mineshige et al. 2000; see also Watarai et al. 2001) 


^Bol 


1 + In 



(for > ^# c ), 


(13) 


where = 2 ~ 3 and ~ 20. The uncertainties on Iq and depend on the details of the 
vertically-averaged equations of slim disks. Calculations of emergent spectra of slim disks show 
that T 5100 oc «^ 2 / 3 M ^ 3 very similar to the SS73 disk, whereas the hydrogen ionizing luminosity 
(Lion = ./'n fjev L e de, where e is photon energy) increases with much slower than L 510 o (Wang 
et al. 1999; Mineshige et al. 2000; Watarai et al. 2001; Wang et al. 2014b) and, for small to 
intermediate size BHs, Li on ~ Lb 0 i for large ./#, showing a characteristic saturation of Li on with 
. A canonic spectral energy distribution, L e oc e _1 exp(—e/eo), is expected from the photon 
trapping part of slim disks, where eo is the cutoff energy determined by the maximum temperature 
of the disk. Such a SED follows from a temperature distribution of T e g oc i?” 1 / 2 , where R is the 
radius. This is significantly flatter than that in the SS73 disk (Wang & Zhou 1999; Mineshige et 
al. 2000). 

The calculations of the structure and emitted spectrum of slim accretion disks are very chal¬ 
lenging because of the complicated coupling between structure and radiation transfer. The classical 
model (Abramowicz et al. 1988) uses the vertically-averaged equations (their validity regimes were 
discussed by Gu & Lu 2007 and Cao & Gu 2015), which probably underestimate photon trapping 
effects (Ohsuga et al. 2002). The model neglects the time delay between energy generation around 
mid-plane and photon escaping from the surface of the disks. Therefore, both the saturated lumi¬ 
nosity and the transition accretion rates are not well determined by theory. Numerical simulations 
of slim disks show different, perhaps more realistic properties, but the agreement between 2D fully 
general relativistic (Sadowski et al. 2014) and 3D magnetohydrodynamic calculations (Jiang et al. 
2014) is not very good. In particular, in 3D simulations, the inclusion of vertical transportation of 
radiation fluxes driven by magnetic buoyancy significantly increases the output power of radiation 
from slim disks. However, such a process could be greatly suppressed by the fast radial motion of 
the accretion flows so that the radiation remains at a level of low efficiency as the classical model 
(Abramowicz et al. 2013). The detailed discussion of these theoretical issues is beyond the scope of 
this paper. Here we focus on the attempt to use our new sample in order to qualitatively compare 
with the saturated luminosity and the transition accretion rates. Quantitative comparison with 
models of slim disks is given in a forthcoming paper. 

The best measured value of Kj a t (Equation 12) results in two expressions for the critical BLR 

size, 

r H/3 = 3 . 3 ^q 5 x 10 4 m 9 ' 29 ; or i? cr it(H/3) = 18. 8+2 § m Y ‘ 29 ltd, (14) 

which are independent of the accretion rates for AGNs with ^ ^ c . In the context of pho¬ 
toionization, the critical BLR size corresponds to a saturation of the ionizing luminosity, and/or 
the bolometric luminosity. In particular, the small scatter of Y^ at indicates that the scatter of the 
saturated luminosity is small in the current SEAMBH sample. Obviously, the Y —saturation qual¬ 
itatively supports the saturated luminosity shown by Equation (13), but Equation (14) provides a 
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quantitative way to test theoretical models and numerical simulations of slim disks. The ultimate 
test for this finding will be performed once there are enough objects with the same BH mass with 
directly measured R H/3 . 

While the uncertainty on the preferred value of the transition accretion rate is large, because of 
the small number of objects, it is clear that the value obtained from fitting our sample is considerably 
lower than = 50 preferred by Mineshige et al. (2000), but covers = 20 within the error bar 
given by Watarai et al. (2001). It has been suggested that the classical slim disk models may have 
underestimated the effects of photon trapping. In fact, < 10 is predicted by model A of Ohsuga 
et al. (2002; their Figure 1). In this model, it is assumed that most of dissipated gravitational 
energy is released in the mid-plane and the vertical structure of gas density is presumed; however, 
this assumption is not self-consistent. We also note that the minimum accretion rates of -M m \ n = 3 
used by us to define SEAMBHs is outside the range found here. However, there must be a smooth 
transition from thin to slim disks and a relatively large range in ^ from linear dependence on 
^0.35 saturation. 

In summary, the Y —saturation measured from RM experiments provides a new tool to test 
the physics of accretion onto black holes. We find that the current data support both sides of disk 
paradigm, from thin disks with the expected dependence on mass and accretion rate, to slim disks 
with their expected saturation. 


5.2.4■ The Maximum Size of the BLR 


We can also calculate the maximum expected dimensionless radius, r H/3 , using pure observa¬ 
tional considerations. For this we use Equation (3) and the definition of r H/3 from Equation (9). 
This gives, r H/3 = i? H/3 /R g = f~^ n (c/Vfwhm) • This allows us to use the smallest observed H f3 line 
width to derive a mass-independent maximum dimensionless size for all RM sources, 


r 


max 

H/3 



9 X !0 4 


(15) 


where = bmm/lO'^kni s 1 is the minimum observed FWHM(H/3). 

Our sample includes one source, Mrk493, with FWHM(H/3) ~ 780km s -1 . This width is, 
arguably, the smallest FWHM observed so far (see e.g., Zhou et al. 2006; Hu et al. 2008b). 
This translates to ~ 1.5 x 10 5 . As explained, this maximum gravitationally scaled size is 
independent of the BH mass. For systems powered by thin accretion disks we use the definition 
of ?’ cr it = 2.7 x 10 4 mf 0 ' 29 (by inserting ^ = 3 into Equation 9) to find the BH mass requires to 
get such a scaled radius, 7777 ~ 0.003. The data in Table 6 shows that for Mrk493 777,7 — 0.14 and 
^ ~ 76, which indicates that this AGN is not powered by a thin disk. It is therefore not surprising 
that the two mass estimates are not in agreement. 
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5.3. Slim Disks and the Size of the Torus 

The emissivity-weighted radius of the BLR measured by RM experiments is directly related 
to dust in the outskirts of the BLR, probably in a toroidal-shaped structure (the “central torus”). 
As argued in Barvainis (1987) and Netzer & Laor (1993), the innermost boundary of the torus 
is set by the dust sublimation temperature which introduces a simple Rdust oc L Bol relationship. 
Recent dust RM experiments (Suganuma et al. 2006; Koshida et al. 2014) provided i?dust hi 
17 sources by following the response of the dust emission in the JL-band to the variable K-band 
continuum. The results are in very close agreement with the estimated graphite-grain sublimation 
radius assuming a mean grain size of about 0.05/mi. Mor & Netzer (2012) estimated this radius 
to be, Rdust = 0.5L^ 2 T, 2 8 f J0 pc, where L 46 = L B oi/10 46 erg s _1 and Ti 800 = T sub /1800K is the 
sublimation temperature. 

A comparison of the Koshida et al. (2014) observations and the ones presented here for ./# < 3 
sources shows that, on average, i?dust/7? H/3 — 4. For ./# > 3 objects, we find Rdust/7? H(9 — 7 (based 
only on four objects). Since the numbers are too small to derive a meaningful correlations for the 
two accretion rate sub-groups, we focus on the comparison of the measured BLR and torus sizes. 

Figure 7 compares the dust and H/3 radii as a function of ./# in two different ways. The left 
panel (7a) shows the dimensionless radii of the BLR and the torus inner walls for 14 objects in the 
Koshida et al. (2014) sample for which we have measured H/3 lags. The right panel (7b) shows the 
ratio of the two radii for all these objects. The surprising result is the unexpectedly large scatter in 
this ratio. While H (3 and the hot dust follow the continuum luminosity in basically an identical way 
(size proportional to L 4 / 2 ), and differ only in their normalization, this does not result in similar 
ratios between the two distances for all sources. In fact, the ratios we measure range from about 
1.4 for Mrk590 to about 10 for Mrk335 and NGC4593. The median ratio, shown in the diagram, 
is about 4. The large scatter and the very small number of sources with measured dust lags prevent 
us from obtaining a meaningful relationship between i?dust and 

There are several interesting ideas about the location of the innermost walls of the torus around 
slim accretion disks (e.g., Kawakatu & Ohsuga 2011; Kawaguchi 2013). In particular, the large 
attenuation of the slim disk radiation at high polar angles (Wang et al. 2014b) would predict 
a small i?dust at these locations, perhaps as small as the self-gravity radius of the disk. This is 
not evident in the present data which suggests that the geometry of the inner torus walls is more 
complicated than the simplistic model suggested in various papers. 


5.4. SEAMBHs and Cosmology 

The new results presented here suggest that the i? H/3 — L 5100 relationship for AGNs requires 
an additional parameter. In particular, the use of the earlier correlations that are based on low 
accretion rate sources, tend to over-estimate the BLR size of high—./# sources and hence their BH 
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mass. The shorter H f3 lags found here indicate that these objects have smaller BLR size than the 
ones obtained by applying the old R H6 — T 5100 relationships. This implies smaller BH mass and, 
therefore, higher Eddington ratios. Recent work by Kelly & Shen (2013), Nobuta et al. (2012), 
and Netzer & Trakhtenbrot (2014) all show the large fraction of high Aboi/TecM sources at high 
redshifts. This number is likely to increase given the necessary correction to the BH mass found 
here. 

There are several ways of using AGNs for cosmology. These include using the measured 
diameters of various BLRs (Elvis & Karovska 2002), the empirical R Hl3 — T 5100 relation (Horne et 
al. 2003; Watson et al. 2011; Czerny et al. 2013), torus inner edge reverberation (Yoshii et al. 
2014; Honig et al. 2014), and RMS X-ray spectra (La Franca et al. 2014). The use of SEAMBHs 
as cosmological distance indicators has been suggested by Wang et al. (2013) and explained in 
great detail in Paper II where we examined the consequences of increasing that leads to an 
almost constant L^ 0 i/M 9 . Marziani & Sulentic (2014) used a similar approach based on a sample 
of highly accreting AGNs identified by their 4D eigenvector 1 (see also King et al. 2014). While we 
used RM-based BH mass measurements (Paper II), they used the standard single-epoch method, 
combined with photoionization modeling, to estimate i? H/3 

The usefulness of the present method for cosmology depends on the number of SEAMBHs 
with accurately measured M, at high redshifts. As shown above, the reduction in the BLR size in 
SEAMBHs compared to the size in low accretion rate sources (used so far to estimate BH mass) 
depends on and can reach a factor of 3 ~ 4 for = 100 (Figure 4a and Equation 7). It 
is reasonable to expect that higher values of are associated with an even greater reduction 
relative to expressions currently used in the literature. At redshifts z = 1 ~ 2, the reduction 
is more significant than the increase in the measured lag due to the cosmological dilution factor 
(1 + z) suggesting that many SEAMBHs at high redshifts, even those with higher luminosity, can 
be monitored, successfully, and their i? H/3 measured, during one observing season, given adequate 
observing facilities. 


6. Conclusions 

We have completed two years of spectroscopic monitoring of 18 AGNs that are high-accretion 
candidates. The second year observations resulted in significant H f3 lag measurements of five AGNs, 
all of which are identified as SEAMBHs. The highest accretion rate reaches > 200 in a few 
objects. The main results of our study can be summarized as follows: 

• For a given L 5100 , the SEAMBHs monitored by us generally have shorter H j3 lags by a factor 
of up to ~ 3 — 4 compared with lower accretion rate AGNs. The reduction gets larger with 
increasing dimensionless accretion rates. 

• We considered several possible explanations for the new results and suggest that they may 
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be related to the strong anisotropy and self-shadowing of the central slim disk radiation. 
We defined a new radius-nrass parameter, Y = m 9 ' 29 r H/3 and demonstrated that it reaches 
saturation when the accretion rate exceeds a critical value of c = 13.8+g 9 { 6 . This transition 
accretion rate is significantly smaller than the one predicted by classical slim disks model, 
but is in better agreement with modified slim disk models. 

• The shorter H (5 lags imply that the number of SEAMBHs in the universe probably has 
been underestimated in earlier works that did not take into account the smaller BLR in 
such sources. It also helps to obtain lag measurements in a shorter time, thus making the 
application of SEAMBHs to cosmology more feasible. 

Our observed sample, combined with SEAMBHs obtained from the literature, is too small to 
put stronger constraints on the radius-mass saturation parameter. Larger SEAMBH samples with a 
large dynamic range in ./#, will provide better understanding of the nature of transition from thin to 
slim disk systems. Better and more accurate BH mass measurements, based on improved numerical 
tools such as the MCMC technique as applied to velocity-resolved RM observations (Pancoast et 
al. 2011; Li et al. 2013), will help too. Finally, we are aware of the limitations of the current theory 
and numerical simulations of slim accretion disks. Present super-Eddington accretion models are 
still missing important physical ingredients and the observations presented here will be compared 
with more detailed calculations once these become available. 
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Table 1. The SEAMBH project: targets and observations 


Object 

CC2000 

<52000 

redshift 

monitoring period 

Aspec 

m 

r'(Lj) 

S/N 

m r /(WO) 

Spec 

Comparison stars 
P.A. 

Note on r H/3 

First phase: SEAMBH2012 sample 

Mrk 335 

00 06 19.5 

+20 12 10 

0.0258 

Oct., 2012 

- Feb., 2013 

91 

50 


21 

00 

O 

174.5° 

Yes 

Mrk 1044 

02 30 05.5 

-08 59 53 

0.0165 

Oct., 2012 

- Feb., 2013 

77 

91 


21 

207".0 

-143.0° 

Yes 

IRAS 04416+1215 

04 44 28.8 

+12 21 12 

0.0889 

Oct., 2012 

- Mar., 2013 

92 

12 


31 

137".9 

-55.0° 

Yes 

Mrk 382 

07 55 25.3 

+39 11 10 

0.0337 

Oct., 2012 

— May., 2013 

123 

25 


3 

198".4 

-24.6° 

Yes 

Mrk 142 

10 25 31.3 

+51 40 35 

0.0449 

Nov., 2012 

- Apr., 2013 

119 

45 


10 

113".l 

155.2° 

Yes 

MCG +06 - 26 - 012 

11 39 13.9 

+33 55 51 

0.0328 

Jan., 2013 

— Jun., 2013 

34 

24 


37 

204".3 

46.1° 

Yes 

Mrk 42 

11 53 41.8 

+46 12 43 

0.0246 

Jan., 2013 

- Apr., 2013 

53 

24 


15 

234".4 

33.8° 

No 

IRAS F12397+3333 

12 42 10.6 

+33 17 03 

0.0435 

Jan., 2013 

— May., 2013 

51 

53 


32 

189".0 

130.0° 

Yes 

Mrk 486 

15 36 38.3 

+54 33 33 

0.0389 

Mar., 2013 

- Jul., 2013 

45 

100 


44 

193".8 

-167.0° 

Yes 

Mrk 493 

15 59 09.6 

+35 01 47 

0.0313 

Apr., 2013 

— Jun., 2013 

27 

40 


45 

155".3 

98.5° 

Yes 

Second phase: SEAMBH2013 sample 

SDSS J075101.42+291419.1 

07 51 01.4 

+29 14 19 

0.1208 

Nov., 2013 

- May., 2014 

38 

50 

50 

40 

133".3 

-41.3° 

Yes 

SDSS J080101.41+184840.7 

08 01 01.4 

+18 48 40 

0.1396 

Nov., 2013 

- Apr., 2014 

34 

50 

50 

50 

118".8 

-98.2° 

Yes 

SDSS J080131.58+354436.4 

08 01 31.6 

+35 44 36 

0.1786 

Nov., 2013 

- Apr., 2014 

31 

50 

50 

30 

100".0 

145.2° 

uncertain 

SDSS J081441.91+212918.5 

08 14 41.9 

+21 29 19 

0.1628 

Nov., 2013 

- May., 2014 

34 

30 

30 

30 

79". 0 

73.9° 

Yes 

SDSS J081456.10+532533.5 

08 14 56.1 

+53 25 34 

0.1197 

Nov., 2013 

- Apr., 2014 

27 

40 


50 

164".5 

-172.9° 

Yes 

SDSS J093922.89+370943.9 

09 39 22.9 

+37 09 44 

0.1859 

Nov., 2013 

- Jun., 2014 

26 

40 

30 

30 

175".1 

-139.0° 

Yes 

SDSS J094422.13+103739.7 

09 44 22.1 

+10 37 40 

0.2410 

Jan., 2014 

- Apr., 2014 

17 

30 

30 

15 

109".3 

170.1° 

No 

SDSS J100055.71+314001.2 

10 00 55.7 

+31 40 01 

0.1948 

Nov., 2013 

- Jun., 2014 

26 

100 

15 

40 

151".2 

116.8° 

No 


Note. — We denote the samples monitored during the 2012-2013 and 2013-2014 observing seasons as SEAMBH2012 and SEAMBH2013, respectively. 
The SEAMBH2012 group is described in Papers I, II and III. No observational points were removed from the SEAMBH2012 sample due to S/N ratios. 
N spec is the numbers of spectroscopic epochs, R * is the angular distance between the object and the comparison star and PA the position angle from the 
AGN to the comparison star. The last column contains notes on the H/3 time lags: “Yes” means significant lag, “No” lags that could not be measured 
and “uncertain” refers to the case of J080131 described in the text. m r /(Lj/WO) are referred to photometry at Lijiang Station of Yunnan Observatory 
and Wise Observatory. We calculated S/N ratios of the SEAMBH2012 and SEAMBH2013 samples from light curves. We removed those points of 
SEAMBH2013 sample in poor observing conditions in light of the lowest S/N ratios listed here. Only a few points with extremely large error bars were 
removed in SEAMBH2012 sample (see details in Papers I and II). Objects marked with “■ • • ” were not observed at the Wise observatory. 






Table 2. Light curves of J075101 and J080101 




J075101 





J080101 



Continuum 

Combined Continuum 


Line 

Continuum 

Combined Continuum 


Line 

JD 

±5100 

JD 

±5100 

JD 

±H/3 

JD 

±5100 

JD 

±5100 

JD 

±H/3 

3.466 

6.219 ±0.169W 

3.466 

6.219 ±0.169 

10.277 

4.312 ±0.055 

13.310 

6.835 ± 0.089L 

13.310 

6.835 ± 0.089 

16.257 

6.696 ± 0.076 

8.588 

5.850 ± 0.181W 

8.588 

5.850 ±0.181 

13.247 

4.413 ±0.074 

16.257 

7.120 ±0.051S 

16.257 

7.120 ±0.051 

21.430 

6.484 ± 0.055 

10.277 

6.325 ± 0.046S 

10.277 

6.325 ± 0.046 

15.248 

4.508 ± 0.091 

21.430 

6.834 ± 0.046S 

21.440 

6.793 ± 0.057 

24.365 

6.294 ± 0.058 

13.247 

6.559 ± 0.075S 

13.256 

6.455 ±0.146 

22.222 

4.019 ±0.053 

21.449 

6.753 ± 0.085L 

24.377 

6.803 ± 0.064 

29.270 

6.330 ± 0.051 

13.265 

6.352 ± 0.103L 

15.248 

7.155 ±0.074 

25.255 

3.699 ± 0.058 

24.365 

6.758 ± 0.076S 

29.283 

6.866 ± 0.060 

33.432 

6.485 ± 0.062 


Note. — JD: Julian dates from 2,456,600; ±5100 and ±h /3 are fluxes at (1 + «)5100A and H/J emission lines in units of 
10“ 16 erg s _1 cm -2 A -1 and 10 -14 erg s 1 cm 2 . L: Photometry from Lijiang. S: Spectrum from Lijiang. W: Photometry from Wise. 
(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its 
form and content.) 




Table 3. Light curves of J080131 and J081441 





J080131 





J081441 




Continuum 

Combined Continuum 


Line 

Continuum 

Combined Continuum 


Line 

JD 

T 5100 

JD 

-F 5100 

JD 

Th/3 

JD 

-F 5100 

JD 

-F 5100 

JD 

Th/3 

3.478 

2.146 ± 0.072W 

3.478 

2.146 ± 0.072 

6.296 

1.301 ±0.022 

3.490 

2.720 ± 0.042W 

3.490 

2.720 ± 0.042 

18.318 

2.590 ± 0.054 

6.296 

2.230 ± 0.020S 

6.311 

2.326 ±0.136 

7.336 

1.318 ±0.021 

8.612 

2.725 ± 0.036W 

8.612 

2.725 ± 0.036 

23.290 

2.856 ± 0.035 

6.325 

2.422 ± 0.029L 

7.351 

2.258 ± 0.040 

9.343 

1.218 ±0.020 

18.318 

2.502 ± 0.064S 

18.318 

2.502 ± 0.064 

30.409 

2.782 ± 0.022 

7.336 

2.287 ± 0.015S 

8.600 

2.146 ± 0.068 

12.401 

1.221 ±0.027 

23.290 

2.599 ± 0.028S 

23.309 

2.626 ± 0.039 

34.369 

2.743 ± 0.028 

7.367 

2.230 ± 0.029L 

9.343 

2.241 ± 0.012 

21.278 

1.175 ±0.024 

23.328 

2.654 ± 0.097L 

30.409 

2.772 ± 0.010 

37.301 

2.728 ± 0.022 


Note. — This table is available in its entirety in a machine-readable form in the online journal, 
regarding its form and content. 


A portion is shown here for guidance 


to 

00 




Table 4. Light curves of J081456 and J093922 




J081456 





J093922 



Continuum 

Combined Continuum 


Line 

Continuum 

Combined Continuum 


Line 

JD 

±5100 

JD 

±5100 

JD 

±H/3 

JD 

±5100 

JD 

±5100 

JD 

±H/3 

6.344 

6.026 ± 0.028S 

6.349 

6.089 ± 0.090 

6.344 

3.507 ± 0.041 

3.589 

2.673 ± 0.086W 

3.589 

2.673 ± 0.086 

17.372 

1.300 ±0.033 

6.354 

6.152 ± 0.031L 

7.380 

6.066 ± 0.030 

7.380 

3.441 ± 0.040 

17.372 

2.668 ±0.014S 

17.372 

2.668 ±0.014 

66.424 

1.137 ±0.029 

7.380 

6.066 ± 0.030S 

10.430 

6.083 ± 0.040 

10.422 

3.297 ± 0.039 

50.596 

2.554 ±0.071W 

50.596 

2.554 ±0.071 

68.304 

1.070 ±0.032 

10.422 

6.054 ± 0.030S 

13.346 

6.017 ±0.046 

19.319 

3.184 ±0.036 

54.414 

2.377 ± 0.078W 

54.414 

2.377 ±0.078 

79.421 

1.029 ±0.038 

10.438 

6.112 ±0.024L 

19.319 

5.974 ±0.051 

22.382 

3.026 ± 0.046 

56.441 

2.537 ± 0.076W 

56.441 

2.537 ±0.076 

81.239 

1.140 ±0.022 


Note. — This table is available in its entirety in a machine-readable form in the online journal, 
regarding its form and content. 


A portion is shown here for guidance 




Table 5. Light curves of J094422 and J100055 




J094422 





J100055 



Continuum 

Combined Continuum 


Line 

Continuum 

Combined Continuum 


Line 

JD 

±5100 

JD 

±5100 

JD 

±H/3 

JD 

±5100 

JD 

±5100 

JD 

±H/3 

3.614 

1.500 ± 0.038W 

3.614 

1.500 ±0.038 

62.387 

1.051 ±0.020 

9.512 

2.921 ± 0.248W 

9.512 

2.921 ±0.248 

18.388 

0.697 ±0.040 

9.476 

1.405 ± 0.063W 

9.476 

1.405 ±0.063 

67.447 

1.043 ±0.027 

18.388 

3.649 ± 0.042S 

18.388 

3.649 ± 0.042 

44.397 

0.823 ± 0.046 

53.437 

1.157 ± 0.035W 

53.437 

1.157 ±0.035 

81.351 

1.030 ±0.029 

43.551 

2.932 ± 0.237W 

43.551 

2.932 ± 0.237 

64.438 

0.747 ±0.030 

56.469 

1.216 ± 0.040W 

56.469 

1.216 ±0.040 

84.394 

0.964 ± 0.024 

44.397 

3.483 ± 0.038S 

44.397 

3.483 ± 0.038 

67.261 

0.798 ± 0.029 

61.528 

1.222 ± 0.040W 

61.528 

1.222 ±0.040 

87.234 

0.961 ±0.018 

45.611 

3.203 ± 0.265W 

45.611 

3.203 ± 0.265 

71.420 

0.802 ± 0.033 


Note. — This table is available in its entirety in a machine-readable form in the online journal, 
regarding its form and content. 


A portion is shown here for guidance 
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Table 6. H/3 Reverberations of the SEAMBH Targets 


Objects 

T H/3 

FWHM 

^line 

log (M./M 0 ) 

log 

log+5100 

lo g l hp 

EW(H/3) 


(days) 

(km s _1 ) 

(km s -1 ) 



(erg s” 1 ) 

(erg s -1 ) 

(A) 

SEAMBH2012 

Mrk 335 

8-7±i;| 

2096 ±170 

1470 ± 50 

6.87i°;i° 

1 oq+0.37 
1.ZO-030 

43.69 + 0.06 

42.03 ± 0.06 

110.5 + 22.3 

Mrk 1044 

10.5^;3 

1178 + 22 

766 + 8 

6-45ig;i3 

1 OO+0.40 
i.zz-o 41 

43.10 + 0.10 

41.39 + 0.09 

101.4 + 31.9 

Mrk 382 

7-5ti;g 

1462 ± 296 

840 ± 37 

6 50" 1 ” 0 ' 19 
D.OU—o 29 

1 1 q+0-69 
1-1O-0.53 

43.12 + 0.08 

41.01 + 0.05 

39.6 + 9.0 

Mrk 142 

7-9tin 

1588 ± 58 

948 ± 12 

6.59i°;°7 

1 «k+0-23 
i.uO—023 

43.56 + 0.06 

41.60 + 0.04 

55.2 + 9.5 

MCG +06-26-012° 

24.0 tit 

1334 ± 80 

785 ± 21 

6-92±g;i| 

—0 34+0.37 

u -o zt _o.45 

42.67 + 0.11 

41.03 + 0.06 

114.6 + 32.5 

IRAS F12397 6 

9-7±5;| 

1802 ± 560 

1150 ± 122 

1 6 7Q + 0 - 27 

D - ' y —0.45 

9 9 ^+ 0 .98 

z.zu_062 

44.23 ± 0.05 

42.26 ± 0.04 

54.2 ± 8.4 

Mrk 486 

23.7±2:? 

1942 ± 67 

1296 ± 23 

7 24 -1-0 ' 12 
' * Z4 _0.06 

0 5 [ 5 + 0 - 20 
U.OO—032 

43.69 + 0.05 

42.12 + 0.04 

135.9 + 20.3 

Mrk 493 


778 ± 12 

513 + 5 

6 14+ 004 

l-oo_0.21 

43.11 + 0.08 

41.35 + 0.05 

87.4 + 18.1 

IRAS 04416° 

13.3ti 3 4 9 

1522 ± 44 

1056 ± 29 

70+0.31 

O- <0-0.06 

2 fi3 + 016 

z.uO-067 

44.47 ± 0.03 

42.51 + 0.02 

55.8 + 4.7 

SEAMBH2013 

SDSS J075101 

33-4^g 5 g 6 

1495 ± 67 

1055 ± 32 

7-i6±g;J5 

1 04 +O .25 
l-04-o 41 

44.12 + 0.05 

42.25 ± 0.03 

68.1 + 8.6 

SDSS J080101 

o q+9-7 
o-o_ 2.7 

1930 ± 18 

1119 ± 3 

6-7818$ 

9 qq+0.39 
Z.OO_072 

44.27 ± 0.03 

42.58 ± 0.02 

105.5 + 8.3 

SDSS J081441 

18.4i,£ 7 

1615 + 22 

1122 ± 11 

6 Q7 + 0 23 
D - y '-0.27 

1 5R + 0,63 
l.OO_ 0 57 

44.01 + 0.07 

42.42 ± 0.03 

132.0 + 23.7 

SDSS J081456 

24.3i 7 6 7 4 

2409 ± 61 

1438 ± 32 

7 44 + 0 - 12 
' *^- 0.49 

0 5Q + 103 
u.oy_ 0 30 

43.99 + 0.04 

42.15 + 0.03 

74.4 ± 7.6 

SDSS J093922 

n-9tS:s 

1209 ± 16 

835 ± 11 

6 53 "*”°' 07 

u.oo_033 

2 e ) 4 + 0,71 
z -°^-0.20 

44.07 ± 0.04 

42.09 ± 0.04 

53.0 + 6.7 


Note. - 1: H/3 lags of J080131 are not listed in this table, but given in the second paragraph of Section 3 in the 

main text. 2: All SEAMBH2012 measurements are taken from Paper III, but 5100A fluxes are from I and II. “MCG 
+06—26—012 was selected as a super-Eddington candidate but later was identified to be a sub-Eddington accretor 
{.M = 0.46). 6 For IRAS F12397, we use the fluxes of the case with local absorption correction (see the details in 
Paper I and III). c The time lag of IRAS 04416 cannot be obtained significantly using the integration method in 
Papers I and II, but has been detected by the fitting procedures in Paper III. 
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Table 7. Reverberation mapping AGNs and results 


Objects 

r H/3 

FWHM 

log (M./Mq) 

log^f 

log A5100 

lo g L H p 

EW(H/3) 

Ref. 


(days) 

(km s -1 ) 



(erg s _1 ) 

(erg s _1 ) 

(A) 


Mrk 335 

8.7±i;® 

2096 ±170 

6.87l?;“ 

1 1 7+0.31 
±-± ' —0.30 

43.69 + 0.06 

42.03 ± 0.06 

110.5 ± 22.3 

1,2,3 


16-8±|;1 

1792 ± 3 

7-021?;“ 

1 OO+0.30 

1 .^ 8 _o 29 

43.76 ± 0.06 

42.13 + 0.06 

119.7 + 23.3 

4,5,6,7 


12 5+ 6 - 6 

1679 ± 2 

6-84l°;“ 

1 39" 1 ” 0 ’ 30 
± 0 -0.29 

43.84 ± 0.06 

42.18 + 0.06 

111.2 ± 21.1 

4,5,6,7 


14 3+ 0,7 

1724 ± 236 

6 92 +0,11 

°* yz -0.14 

1 9 c+0.30 
1ZO -0.29 

43.74 ± 0.06 

41.99 + 0.07 

89.5 ± 19.5 

4,8“ 


i4.o ±t:t 

1707± 79 

6 93 “*" 010 

o.w«3_o 11 

1 97 +O.I 8 

‘ — 0.17 

43.76 ±0.07 

42.09 ±0.09 

108.2 ± 16.8 


PG 0026+129 

lll.ot^.i 

2544 ± 56 

0 ic+0.09 

0 6 C )+ 0 - 28 
u. 00-0 20 

44.97 + 0.02 

42.93 ± 0.04 

46.2 + 4.7 

4,5, 6 ,9 

PG 0052+251 

89.81**;* 

5008 ± 73 

8.64l?;“ 

—0 59 ~*~°' 31 
u,oy -0.25 

44.81 + 0.03 

43.13 + 0.05 

107.4 + 14.8 

4,5, 6 ,9 

Fairall9 

17 4"*" 3,2 

1 ' '^—4.3 

5999 ± 66 

8.091?;?? 

—0 7l + 0 - 31 
u - 1 -0.21 

43.98 ± 0.04 

42.67 + 0.04 

249.8 ± 32.0 

4,5,6,10 

Mrk 590 

20.7I2;* 

2788 ± 29 

7.501?;?? 

—0 22" 1 " 0 ' 24 
U-ZZ -0.25 

43.59 ± 0.06 

41.92 + 0.06 

107.0 + 22.0 

4,5, 6 ,7 


14 0 " 1 " 8 ' 5 
8.8 

3729 ± 426 

7 5 «+ 0 - 22 

<-Oo _ 0 48 

—0 qi + 0 - 28 
0.30 

43.14 + 0.09 

41.58 + 0.16 

142.8 ± 62.2 

4,5, 6 ,7 


29 2 ~*~ 4 ' 9 

2744 ± 79 

7 63”*” 0 ' 07 

1 -VO_Q 09 

_ 0 54+0-25 

U-O^_o.26 

43.38 ± 0.07 

41.75 + 0.07 

119.7 + 28.1 

4,5, 6 ,7 


28 . 811 ® 

2500 ± 42 

7 cc+0.05 
' - oo -0.07 

—0 13" 1 " 0 ' 24 
u,±o -0.25 

43.65 ± 0.06 

41.92 + 0.07 

94.9 + 20.7 

4,5, 6 ,7 


25.6l®| 

2716 ± 202 

7 rr+0.07 
' * DD —0.08 

—0 41+ 0 * 36 
u, ^ t —0.36 

43.50 ±0.21 

41.85 ±0.12 

108.6 ± 20.2 


Mrk 1044 

10.51®;? 

1178 ± 22 

6 45 +0 ' 12 
o- z ±o _ 0 13 

1 99+0.40 
i.zz-o 41 

43.10 + 0.10 

41.39 + 0.09 

101.4 + 31.9 

1,2,3 

3C 120 

38.il?i;i 

2327 ± 48 

7 61" 1 ” 0 ' 19 

1 - d± -0.22 

0 03+ 0 - 37 

O.UO -0 37 

44.07 + 0.05 

42.37 + 0.06 

100.9 ± 18.3 

4,5, 6 ,7 


25.91?;? 

3529 ± 176 

' -° U -0.06 

—0 17+ 0 - 37 
u ± ' -0.37 

43.94 ± 0.05 

42.36 ± 0.05 

135.9 ± 22.3 

4,8“ 


26 - 21 *- 

2472 ± 729 

7 79+ 015 

• * * y — 0.15 

1 

O 

b 

1 + 
0 0 

w co 
0 O 

44.00 ±0.10 

42.36 ±0.04 

118.8 ±28.9 


IRAS 04416 

13.3l? 3 4 9 

1522 ± 44 

70+0.31 

O. ( o_o .06 

2 53 + 0 - 16 

z.uo-o 67 

44.47 + 0.03 

42.51 + 0.02 

55.8 + 4.7 

1,2,3 

Ark 120 

47 1 + 8.3 

6042 ± 35 

8.531?;?? 

1 4&+0 ,2 4 
t.4o_o 23 

43.98 ± 0.06 

42.60 ± 0.05 

211.5 ± 37.5 

4,5, 6 ,7 


37 i + 4 -8 

6246 ± 78 

8.451?;?? 

— 2 01 + 0,27 
z,ul -0.27 

43.63 ± 0.08 

42.43 ± 0.07 

321.1 ± 77.6 

4,5, 6 ,7 


39.5I?;® 

6077 ± 147 

8 47+°’ 07 

0 ' —0.08 

-1 70+°- 41 
A *' u — 0.41 

43.87 ± 0.25 

42.54 ±0.13 

244.8 ±80.3 


Mrk 79 

9-ol?;i 

5056 ± 85 

7 R 5+ 0 - 28 

1 .uo_ 0.88 

-0 75 + 0 - 41 
u - /o -0.34 

43.63 ± 0.07 

41.89 + 0.07 

92.4 + 21.0 

4,5, 6 ,7 


ie.it®;® 

4760 ± 31 

7 oc+0.15 

/ •oo _ 0< 23 

—0 5Q ~*~ 0,41 
u - oy -0.34 

43.74 ± 0.07 

41.92 + 0.08 

78.1 ± 18.3 

4,5, 6 ,7 


16-0l|;| 

4766 ± 71 

7 or+0.15 

1 .Od _0 20 

-0.70l?;“ 

43.66 ± 0.07 

41.89 + 0.07 

86.0 + 19.3 

4,5, 6 ,7 


15-61*- 

4793 ± 145 

7 84+° 12 

1 —0.16 

- 0 . 681 ?;*? 

43.68 ±0.07 

41.90 ±0.05 

85.4 ± 13.3 


SDSS J075101 

33.4l ? 5 6 6 

1495 ± 67 

7 16”*~ 0 ' 17 

<-lO_ 0 .09 

1 Q4+0.25 
l-o4_o 41 

44.12 + 0.05 

42.25 ± 0.03 

68.1 ± 8.6 

11 

Mrk 382 

7 c+2.9 

1462 ± 296 

6.501?;?? 

1 1 o+0.69 
1.1»_0.53 

43.12 + 0.08 

41.01 + 0.05 

39.6 ± 9.0 

1,2,3 

SDSS J080101 

0 0 + 9.7 

o.O_27 

1930 ± 18 

6.781?;?? 

9 qq + 0.39 
Z.OO-Q 72 

44.27 + 0.03 

42.58 ± 0.02 

105.5 + 8.3 

11 
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Table 7—Continued 


Objects 

r H/3 

FWHM 

log (M./M 0 ) 

logoff 

log A5100 

lo g L H/3 

EW(H/3) 

Ref. 


(days) 

(km s -1 ) 



(erg s -1 ) 

(erg s” 1 ) 

(A) 


PG 0804+761 

146.9li*;® 

3053 ± 38 

O 30 + 0.05 
o-4o_ 0 06 

0.00+(i® 

44.91 + 0.02 

43.29 ± 0.03 

122.5 ± 10.3 

4,5,6,9 

SDSS J081441 

18.4t& 7 

1615 ± 22 

6 Q7+ 0 - 23 

D,y '-0.27 

1-56+1? 

44.01 + 0.07 

42.42 ± 0.03 

132.0 + 23.7 

11 

SDSS J081456 

24-3 + 7 16 7 a 

2409 ± 61 

7 44+ 012 
* ‘”—0.49 

0 59“*" 1,03 
u,oy -0.30 

43.99 ± 0.04 

42.15 + 0.03 

74.4 ± 7.6 

11 

PG 0844+349 

90 q+13.7 

oz.o_ 13 4 

2694 ± 58 

'~7 £?/? + 0.15 

7.66_o.23 

0-50+S 

44.22 + 0.07 

42.56 + 0.05 

111.2 ± 22.1 

5,6,9,12 

Mrk 110 

24.3+|;| 

1543 ± 5 

7 Q5+ 0 - 09 
' - uo -0.18 

0-81+1 

43.68 ± 0.04 

42.12 + 0.05 

139.6 + 20.4 

+5,6,7 


20-4+g°3 5 

1658 ± 3 

7 04 +0 - 18 
0.16 

0 92 +0 ' 34 
u,yz —0.32 

43.75 + 0.04 

42.02 + 0.05 

94.8+ 14.7 

+5,6,7 


qq q+14.9 
oo -°-10.0 

1600 ± 39 

7 qq+ 0.16 
' ,zz —0.16 

0 '58”*”°' 35 
u - oo -0.33 

43.53 + 0.05 

41.97 + 0.04 

139.4 + 20.5 

+5,6,7 


25 . 

1634 ± 83 

7 io+° 13 

1 —0.14 

0 77+ 0 - 26 

U ‘ ' * —0.25 

43.66 + 0.12 

42.03 + 0.08 

123.8 + 29.1 


SDSS J093922 


1209 ± 16 

6 '53+ 0 - 07 
O-OO-Q 33 

2 54 ”*”°' 71 
z, ° — 0.20 

44.07 + 0.04 

42.09 + 0.04 

53.0 + 6.7 

11 

PG 0953+414 

150. 1 + 22(6 

3071±27 

8 44+ 006 
0.07 

0 39 “*“°* 16 

u.oy_ 0 14 

45.19 + 0.01 

43.29 + 0.04 

64.7 + 5.9 

4,5, 6 ,9 

NGC 3227 

q 0+O.8 

4112 + 206 

7 09 +0 09 
' ,u - 0.12 

— 1 34+ 038 
0.36 

42.24 + 0.11 

40.38 + 0.10 

71.0 + 23.6 

4,13“ 

Mrk 142 

7 - 9 ± i :? 

1588 ± 58 

6 59+ 0 - 07 

o.oy_ 0 07 

1 an+°- 85 
j..»u_o 86 

43.56 + 0.06 

41.60 + 0.04 

55.2 ± 9.5 

1,2,3 


9 7+0.7 
' -0.8 

1462 ± 2 

6 06 +0 10 
°- uq -0.16 

1 96”*~°' 82 
±.»u_o 82 

43.61 + 0.04 

41.66 + 0.05 

57.6 ± 8.6 

4,14 


6 4+ 73 

1462 ± 86 

6 47+ 0 - 38 

0 ' —0.38 

1 93“*"°‘ 59 

43.59 + 0.04 

41.62 + 0.06 

56.6 ± 6.6 


NGC 3516 

11 7+!0 

5384 ± 269 

7 09+0.05 

' -o^_o.08 

-1-97+1 

42.79 + 0.20 

41.06 + 0.18 

94.7 + 59.2 

4,13“ 

SBS 1116+583A 

9 q+0.6 

z -°-0.5 

3668 ± 186 

6.78+“ 

—0.87+®i 

42.14 + 0.23 

40.70 + 0.07 

186.8+ 104.1 

4,14 

Arp 151 

4.0t°;® 

3098 ± 69 

/? 07+O.O5 
O-O' -0.08 

—0 44+ 030 
u *”— 0.28 

42.55 + 0.10 

40.95 + 0.11 

130.0 + 44.4 

4,14 

NGC 3783 

10 . 2 + 2(3 

3770 ± 68 

7 . 45 +“ 

1 tro + 0.45 
— i.OO_o 59 

42.56 + 0.18 

41.01 + 0.18 

144.0 + 83.7 

4,5,6,15 

MCG +06-26-012 

24.0l®+ 

1334 ± 80 

6-92+“ 

—0 34 + 0 - 37 
u -o4_o 45 

42.67 + 0.11 

41.03 + 0.06 

114.6 + 32.5 

1,2,3 

Mrk 1310 

q 7+O.6 

2409 ± 24 

6 62"*” 0 07 
0.08 

—0 31 “*~ 0,35 
u *° — 0.39 

42.29 + 0.14 

40.56 + 0.10 

94.3 ± 38.2 

4,14 

NGC 4051 

+9io,5 

851 ± 277 

5 42”*” 0,23 

0,4 —0.53 

1 59 “*" 1,29 
0.84 

41.96 + 0.19 

40.19 + 0.18 

86.8 + 51.7 

4,13“ 

NGC 4151 

6 - 6 tS;J 

6371 ± 150 

7 7Q+0.07 
' • ‘ z — 0.06 

—2.81 + 57 

42.09 + 0.21 

40.56 ± 0.20 

150.8 ± 100.6 

4,5,6,16 

PG 1211+143 

93.812(6 

2012 ± 37 

7 07+O.II 
* -0.26 

0-84+1 

44.73 ± 0.08 

43.02 + 0.06 

100.2 ± 22.9 

5,6,9,12 

Mrk 202 

3-0+I 

1471 ± 18 

6 11+ 0 - 20 
D - 11 - 0.20 

0 66”*”°' 59 

42.26 + 0.14 

40.40 + 0.09 

70.6 ± 27.5 

4,14 

NGC 4253 

6 2 +1 ' 6 
o.z_i 2 

1609 ± 39 

0 49+ 010 
O-^-O.lO 

0 36“*" 0,36 
U.OO-o 42 

42.57 + 0.12 

40.77 + 0.12 

81.1 + 31.6 

4,14 

PG 1226+023 

306.81®®;® 

3509 ± 36 

Q Q7+0.09 
8 - 8 '-0.15 

0 70+°- 33 
u - ‘ u - 0.20 

45.96 + 0.02 

44.13 + 0.04 

74.6 ± 6.9 

4,5, 6 ,9 

PG 1229+204 

37-81*5;® 

3828 ± 54 

8 . 03 + 0(23 

— 1 03”*”°' 52 

l- u O_o.55 

43.70 + 0.05 

42.31 + 0.06 

209.7 + 38.3 

4,5, 6 ,9 

NGC 4593 

q 7+O.8 
' —0.8 

5143 ± 16 

7 qq+0.08 

'•^o_0.10 

—0 73 + 0 - 41 
U -' °-0.52 

42.87 + 0.18 

41.17 + 0.18 

101.6 + 59.0 

4,6,17 
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Table 7—Continued 


Objects 

T H/3 

FWHM 

log (M./M q ) 

log^f 

log £5100 

lo g L H/3 

EW(H/3) 

Ref. 


(days) 

(km s -1 ) 



(erg s” 1 ) 

(erg s” 1 ) 

(A) 



4 q+1-3 

^-o.s 

4395 ± 362 

7 91 +0.13 
' >z± -0.12 

— 1 47+0-41 

1,4 ‘ -0.52 

42.38 + 0.18 

40.73 ± 0.18 

115.4 + 67.6 

18® 


4 0 +°- 8 
—0.7 

5142 ± 572 

7 26^ -0 ' 09 
* —0.09 

—1.10 +0,6 ° 
x —0.64 

42.62 + 0.37 

40.95 + 0.33 

108.3 + 45.7 


IRAS F12397 

9-7 til 

1802 ± 560 

6 79+ 0 - 27 

D - ‘ y -0.45 

9 9f?+0.98 

z,.z,u_ 0 62 

44.23 ± 0.05 

42.26 ± 0.04 

54.2 ± 8.4 

1,2,3 C 

NGC 4748 

5.5 i|;| 

1947 ± 66 

a -i +0.11 

5.51-0.23 

o-io±g;2i 

42.56 + 0.12 

40.98 + 0.10 

136.8 + 50.1 

4,14 

PG 1307+085 

105.6i^ ; o 

5059 ± 133 

q 79+0.13 

1 z -0.26 

-0.68i°;®| 

44.85 ± 0.02 

43.13 ± 0.06 

98.4 + 15.1 

4,5,6,9 

NGC 5273 

o 9+1-2 
z ' z -1.6 

5688 ± 163 

7 14+ 0 - 19 
' 0.56 

-2.501J;®? 

41.54 + 0.16 

39.74 + 0.11 

82.2 + 37.1 

19 

Mrk 279 

16 7+ 3 - 9 
±D -' -3.9 

5354 ± 32 

7 Q7+0.09 
' * -0.12 

-0.89±g;|g 

43.71 + 0.07 

42.12 ± 0.06 

132.2 ± 28.7 

4,5,6,20 

PG 1411+442 

124.3tgJ;° 

2801 ± 43 

O 90 + 0.17 

o-^o_0.30 

r» 90+O.63 
u.z,o_o 38 

44.56 ± 0.02 

42.85 ± 0.03 

99.7 + 8.2 

4,5,6,9 

NGC 5548 

19-7111 

4674 ± 63 

7 Q9+0.03 

1 -»^_o.04 

— 1 62 +0 ' 46 
l.DZ-0 49 

43.39 + 0.10 

41.79 ± 0.10 

128.1 + 40.3 

4,5,6,21 


18-6ll;J 

5418 ± 107 

0 rjo+0.05 
o-Uo_ 0 06 

— 1 qq+°- 47 

J-.»»_o.5i 

43.14 + 0.11 

41.61 ± 0.13 

151.3 ± 57.9 

4,5,6,21 


15-911? 

5236 ± 87 

7 93+ 007 
' ' y °—0.08 

-L.oo-o 49 

43.35 ± 0.09 

41.72 ± 0.10 

119.7 + 37.9 

4,5,6,21 


n-oll? 

5986 ± 95 

7 89"*”°' 07 

1 -o»_o.09 

—2 10 + °' 47 
z,iU -0.52 

43.07 + 0.11 

41.52 ± 0.17 

144.5 ± 66.7 

4,5,6,21 


!3.0ll® 

5930 ± 42 

7 95+ 0,05 
' • yo — 0.05 

1 79+O.46 
a *'^-0.49 

43.32 + 0.10 

41.75 ± 0.09 

135.9 + 41.2 

4,5,6,21 


13 4+ 3 - 8 
±0 ^-4.3 

7378 ± 39 

8 1 '5+ 0 - 11 
°- lo -0.17 

-1 64 +0 - 45 
1 -°^-0.49 

43.38 ± 0.09 

41.73 ± 0.10 

114.4 + 37.0 

4,5,6,21 


9-1 rr+2.6 
'-2.6 

6946 ± 79 

0 qi +0.05 
0.06 

— 1 43+0- 45 
1.40-0.48 

43.52 ± 0.09 

41.82 ± 0.09 

102.4 ± 30.2 

4,5,6,21 


WAtli 

6623 ± 93 

8 1 5+°- 93 
°- lo -0.03 

— 1 'SB - * -0 ' 45 

1.0U-Q48 

43.43 ± 0.09 

41.75 ± 0.10 

106.3 ± 33.3 

4,5,6,21 


17-511? 

6298 ± 65 

q 1 q+0.05 
o - lo -0.04 

1 qc+0.46 
i‘°°—0.49 

43.24 + 0.10 

41.72 ± 0.10 

153.5 ± 50.7 

4,5,6,21 


26.5111 

6177 ± 36 

8-30l?;°I 

1 qq+0.45 
i‘°°—0.48 

43.59 ± 0.09 

41.87 + 0.10 

98.1 + 30.0 

4,5,6,21 


24.8111 

6247 ± 57 

Q 9Q + 0.05 

o-^o_0.06 

_1 4 c+0.45 
1.40_o.48 

43.51 + 0.09 

41.83 ± 0.09 

106.6 ± 31.5 

4,5,6,21 


6.511:? 

6240 ± 77 

7.69t?:?? 

—2 04"*" 0,47 
Z ‘ U 4-0.51 

43.11 + 0.11 

41.64 + 0.13 

172.8 + 66.1 

4,5,6,21 


14.31® 3 

6478 ± 108 

8.07t?;?? 

—2 03 + °' 47 
z.uo_ 0 51 

43.11 + 0.11 

41.55 ± 0.14 

139.8 ± 55.7 

4,5,6,21 


6.3il:i 

6396 ± 167 

7 70"*" 0,15 

1 • /u -0.20 

9 97+O.49 
' —0.55 

42.96 + 0.13 

41.12 ± 0.10 

74.4 ± 27.3 

4,22 


4-2l ?;? 

12771 ± 71 

q 1 9+O.O8 

o-^-o.ie 

—2 19 +0 ' 47 
^•iy_0.51 

43.01 + 0.11 

41.33 ± 0.10 

105.3 ± 35.4 

4,14 


12.4±2;| 

11481 ± 574 

8.50t?;?? 

9 91 +0.48 
-0.53 

42.99 + 0.11 

41.27 + 0.10 

96.6 + 34.4 

4,13° 


17.611;? 

7241 ± 2200 

8.io±g;iS 

1 sn - ^ -9,33 

l.OU — o 33 

43.29 + 0.20 

41.64 + 0.23 

117.2 ± 25.9 


PG 1426+015 

95.011?;? 

7113 ± 160 

8.97t?;?? 

— 1 51+0. 47 
i‘Oi— 0.28 

44.63 ± 0.02 

42.83 ± 0.04 

80.1 ± 9.2 

4,5,6,9 

Mrk 817 

io.ol?;? 

4711 + 49 

7 Q9+O.O8 

' 0.09 

-0.8l±°;f® 

43.79 ± 0.05 

42.07 + 0.05 

98.0 + 16.5 

4,5,6,7 


15.311;? 

5237 ± 67 

7 9i+ 0 - 09 
<-yJ-_o.il 

-0.98tg;“ 

43.67 + 0.05 

42.00 + 0.06 

108.5 ± 20.4 

4,5,6,7 
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Table 7—Continued 


Objects 

T H/3 

FWHM 

log (M./Mq) 

log^f 

log L 5100 

l0 g L H/3 

EW(H/3) 

Ref. 


(days) 

(km s -1 ) 



(erg s _1 ) 

(erg s _1 ) 

(A) 



33.6t®j 

4767 ± 72 

0 1 7+O.O8 

°- i7 -o.n 

-0.981+1! 

43.67 + 0.05 

41.92 + 0.05 

91.1 ± 15.3 

4,5,6,7 


w.oty 

5627 ±30 

7 Q4+0- 09 
' - y ^-0.12 

fl 70+0.35 

—u. IO_ 0 35 

43.84 + 0.05 

41.77 + 0.05 

43.2 ± 7.1 

4,13 


19 . 9 +®;® 

5348 ± 536 

7 99+° 14 

q q 
d d 
+ 1 

b- 

oq 

d 

1 

43.74 ±0.09 

41.93 ±0.14 

78.5 ±34.3 


Mrk 1511 

r --r-1-0.9 
' -0.8 

4171±137 

7 oq+0-07 
' - zy -0.07 

—0 34+O 24 
u -' 3 ^-0.24 

43.16 + 0.06 

41.52 + 0.06 

115.5 ± 23.1 

18 b 

Mrk 290 

0 -7+1.2 

°- 7 -i.o 

4543 ± 227 

7 cc+0.07 
' -°°-0.07 

—0 «5+0- 23 
u>oo -0.23 

43.17 + 0.06 

41.64 + 0.06 

153.0 + 29.0 

4,13° 

Mrk 486 

23 . 7 ^ 2 ;? 

1942 ± 67 

7 24+0- 12 
<-^_0.06 

0 55+0- 20 
U.OO-032 

43.69 ± 0.05 

42.12 + 0.04 

135.9 ± 20.3 

1,2,3 

Mrk 493 


778 ± 12 

6 14+ 004 
D - i4 -0.ll 

1 SS+0-03 

t-oo_ 0> 21 

43.11 + 0.08 

41.35 + 0.05 

87.4 + 18.1 

1,2,3 

PG 1613+658 

40-i±i®;° 

9074 ± 103 

0 01 +0.14 

—0 Q7+0- 45 
u - y ‘ -0.31 

44.77 + 0.02 

43.00 ± 0.03 

86.7 + 7.6 

4,5,6,9 

PG 1617+175 

71 c+29.6 

6641 ± 190 

8.791+H 

-i.soll;!! 

44.39 ± 0.02 

42.74 ± 0.05 

114.8 ± 15.1 

4,5,6,9 

PG 1700+518 

251.8 ±H;| 

2252 ± 85 

8.401+1! 

LOSl+H 

45.59 + 0.01 

43.78 ± 0.02 

78.9 + 4.5 

4,5,6,9 

3C 390.3 

23 - 6 -6.7 

12694 ± 13 

Q O7 + O.IO 
8 ' 8 '-0.15 

q qc+0.60 

o.oo-o 65 

43.68 + 0.10 

42.29 ± 0.05 

206.2 ± 50.7 

4,5,6,23 


46-4 ±H 

13211 ±28 

9 20”*” 0 03 

».z,u_o 03 

9 1 O+0-51 
^ -0.51 

44.50 ± 0.03 

42.78 ± 0.04 

97.1 ± 10.0 

4,24 


44-5 

12796 ±361 

9 1«+°* 23 

y — 0.23 

^• b ^_0.96 

44.43 ±0.58 

42.60 ±0.35 

108.8 ±58.8 


KA 1858+4850 

!3.5+23 

1820 ± 79 

6 04+0.07 
0.09 

0 75+0- 25 

U - * °-0.21 

43.43 ± 0.05 

41.89 + 0.04 

146.9 ± 21.1 

25 d 

NGC 6814 

6 6 +0 ' 9 
°-°-0.9 

3323 ± 7 

7 1 g+0.05 
' • iO _0 . 06 

-1 64 +0 - 46 

42.12 + 0.28 

40.50 + 0.28 

121.6 + 112.2 

4,14 

Mrk 509 

79.61+4 

3015 + 2 

0 1 c+0.03 
°- lo -0.03 

—0 52+0- 13 
u ‘ oz -0.14 

44.19 + 0.05 

42.61 + 0.04 

132.7 + 19.1 

4,5, 6 ,7 

PG 2130+099 

9-6+11 

2450 ± 188 

7 Q5+O.O8 

'- UD -o.io 

L69l+1! 

44.20 ± 0.03 

42.65 ± 0.03 

142.1 ± 13.6 

4,8“ 

NGC 7469 e 

10.8111 

4369 ± 6 

7 60"*” 012 

/.OU_ 0 06 

0.9011!? 

43.51 + 0.11 

41.60 + 0.10 

63.0 + 21.0 

4,26 


4 5+0-7 
4 -°- 0.8 

1722 ± 30 

6 42+0- 06 
o-^_0.09 

0-6311!! 

43.32 + 0.12 

41.71 + 0.09 

124.9 + 44.3 

27,5 


fi K+ 5 - 8 

b * —3.0 

4343 ± 2859 

6 92“*"°' 84 

°- y ^_0.84 


43.43 ±0.15 

41.66 ±0.10 

86.9 ±47.1 



Note. — Ref: 1. Paper I; 2. Paper II; 3. Paper III; 4. Bentz et al. (2013); 5. Collin et al. (2006); 6 . Kaspi et al. (2005); 7. 
Peterson et al. (1998); 8 . Grier et al. (2012); 9. Kaspi et al. (2000); 10. Santos-Lleo et al. (1997); 11. This paper; 12. Bentz et 
al. (2009a); 13. Denney et al. (2010); 14. Bentz et al. (2009b); 15. Stirpe et al. (1994); 16. Bentz et al. (2006); 17. Denney et al. 
(2006); 18. Barth et al. (2013); 19. Bentz et al. (2014); 20. Santos-Lleo et al. (2001); 21. Peterson et al. (2002) and references 
therein; 22. Bentz et al. (2007); 23. Dietrich et al. (1998); 24. Dietrich et al. (2012); 25. Pei et al. (2014); 26. Peterson et al. 
(2014); 27. Collier et al. (1998). 

Superscript a means that the literature does not provide FWHMs from mean spectra or H j3 fluxes with narrow-line components 
subtracted. For these, we scanned and digitized the mean spectra from the papers, then calculated those numbers ourselves, b 
means that the AGN continuum is obtained through spectral fitting decomposition (Ref. 18). c for IRAS F12397, we use fluxes 
for the case with local absorption correction (see the details in Papers I and III), d means that the contribution of host galaxy 
is provided in Ref. 25. e the virial products of the two independent campaigns are quite different, yielding a large error bars on 
the mean M« and ^. In Appendix B, NGC 7469 is shown as an SEAMBH from Collier et al. (1988) and Peterson et al. (2014). 
Both campaigns have high quality data. Numbers in boldface are the weighted averages of all the measurements of this object. 
The details of the average scheme are given in the main text. 
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JD - 2456600 (days) Time lag (days) 



20 40 60 80 100 120 140 160 180 -20 0 20 40 60 

JD - 2456600 (days) Time lag (days) 


Fig. 1. — Light curves and cross correlation results. Each object has six panels: (a, b, c) are light curves of the 
inter-calibrated continuum, combined 5100A continuum, and H/3 emissions, respectively; ( d, e, f) are auto correlation 
function (ACF) of the combined continuum, cross correlation function (CCF) of the combined continuum and H/3 
emission and the Monte-Carlo simulations of peaks (red) and centroid (blue) of lags, respectively. In panels d and 
e, the solid lines show the results of ICCF method and the points with error bars are from ZDCF (Z —transformed 
discrete correlation function). F 5100 and F^p are in unit of 10~ 16 erg s -1 cm -2 A -1 and 10~ 14 erg s _1 cm -2 for all 
quasars, respectively. In panel (a) of J075101, the red dots, green squares and blue diamonds are fluxes from Lijiang 
spectroscopy, photometry and Wise photometry, respectively. The colour points in other objects have the same 
meanings with that of J075101. Bars with terminals as systematic errors are plotted in the corners of the panels (see 
Paper I for details). 
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Figure 1 continued. 
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Figure 1 continued. 
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Figure 1 continued. 
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Figure 1 continued. 
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hgj' logEW(H/9) 


Fig. 2. — The distribution of dimensionless accretion rates of the mapped AGNs and quasars. The previous sample 
refers to the sample of all mapped AGNs summarized by Bentz et al. (2013) (the repeated monitored AGNs are 
regarded as individual ones), and adds four recently mapped AGNs: NGC 7469 updated by Peterson et al. (2014), 
KA 1858+4850 (Pei et al. 2014), Mrk 1511 (Barth et al. 2013) and NGC 5273 (Bentz et al. 2014). The dashed line 
is for the previous sample and the solid is for the SEAMBH sample. Our campaign selected those candidate sources 
with extremely high accretion rates > 10). Though some sources monitored previously have ~ 10, most of 
them have lower accretion rates. Right panel shows the distribution of EW(H/3), showing that SEAMBH sample 
tends to have low EW(H/J). 
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Fig. 3. — The R H p —1/5100 and R H/3 — L H/3 relationships (panel a, c, e) and the deviations of R ap from the regression 
relationships (panel b, d, f). The dotted is the fit for the ^ < 3 sample whereas the dashed is the fit for the ^ > 3 
sample in all panels. Panels (a, b) show the comparison of the previous sample and the SEAMBH campaigns. Clearly 
SEAMBHs are below the R H/3 — L 5100 relation, and the scatter of the relation increases with the inclusion of the 
SEAMBH2012/2013 samples (see the intrinsic scatter given by the numbers below Equations 4 and 5). Panel b shows 
ARh /3 = log (R H/3 /R H p tR _ L ), where R nf) R _ L is given by Equation (4b). Panels ( c,d ) show the two samples of > 3 
and < 3 AGNs with averaged RM values for the R H(3 — L 5100 relation and the residual of Af? H(3 , where R H/3 R _ L 
is given by Equation 5b. Panels (e, f) show the R ap — L H/3 and Arelations. 
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Number Number 

Fig. 4. — Panels (a, b) and (c, d) show the Ai? H/3 — relation and A R ap distributions for f? H/3 — L5100 and 
R H p — L n/3 relations, respectively. The dotted lines indicate ^ = 3 in panel a and c. Both panels show that 
A_R h;3 ~ 0 up to ^ ~ 3, beyond which | Af? H/3 1 increases with This implies that, relative to normal AGNs of the 
same 5100A luminosity, SEAMBHs have H/3 lags that become systematically shorter with increasing accretion rate. 
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M Number jtft Number 

Fig. 5. — Black hole mass-scaled BLR and the radius-mass parameter versus dimensionless accretion rate. In each 
panel, there are three plots to show their relation and deviations. Correlations and probabilities of KS tests are 
marked in the plots. Panel a shows that the < 3 sample has a steeper correlation than that given by Equation 
(9). Panel b shows that SEAMBHs deviate from this relation and that Y is independent of Clearly Y behaves 
in a different way in the two accretion rate groups. 
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Fig. 6. — Determination of the saturated —Y and the transition accretion rate. We find Y oc ..#°' 2 ' ±0 ' 04 , supporting 
the SS73 disk model below and saturation beyond the numbers listed inside the panel supporting the idea of slim 
accretion disks in SEAMBHs. The outlier with large error bars is NGC 7469 (see Table 7). 
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Fig. 7. — Black hole mass-scaled BLR (purple and grey points) and dust reverberation radii (black points). Panel 
a shows -Rdust/Rg and R n/3 /Rg, versus ,/M . Panel b shows ratios of Rdust/R H/3 as a function of . Note the large 
scatter in both low and high groups. 
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A. Mean and RMS spectra 


To supplement the material presented in the body of the paper, we provide in Figure 8 the 
mean and RMS spectra of the objects in the SEAMBH2013 sample calculated by the standard 
method (e.g., Peterson et al. 2004). Prior to the calculation of these spectra, we corrected the data 
for the slight wavelength shift caused by the relatively wide aperture used in 2013 (5 W ) using the 
[O m]A5007 line as our wavelength reference. We define 


1 N 


i= 1 


and 


N 


1/2 


s ( a Hv£[ f -w- f w ] 2 


N 


(Al) 


(A2) 


i —1 


where F t (X) is the i —th spectrum and N the total number of spectra obtained during the campaign. 
The line dispersion is 

°Une(^) = A) — ^0 j (A3) 

where Ao = f XP(X)dX/ f P(X)dX, (A 2 ) = f X 2 P(X)dX/ f P(X)dX, and P( A) is the line profile. The 
calculated values of <7ii ne are listed in Table 6. 


B. R h/3 Correlations using the direct method 

The correlations listed in the main body of the paper were calculated using the average scheme 
where each object is represented by one point. The results shown in this section are for the direct 
scheme, where each observing campaign is represented by one point. 

The following equations correspond to Equation (6) 


l°g (^W ltd ) 


' (1.29 ± 0.02) + (0.52 ± 0.03) log L H /j ,42 
< (1.35 ± 0.03) + (0.53 ± 0.03) log L H/ j ,42 
(1.17 ± 0.05) + (0.55 ± 0.06) log L n , 3A2 


(entire sample), (Bla) 

(A < 3), (Bib) 

{A > 3), (Blc) 


with scatters of a m = (0.20,0.15,0.24) for (Bla,Bib,Blc), respectively. For A R H/3 — ^ relation, 
we have 

= (0.49 ± 0.14) - (0.58 ± 0.10) log A (for A > 3), (B2) 

and 

AR Hff = (0.44 ± 0.15) - (0.49 ±0.11) log A (for A > 3), (B3) 

with (Ti n = (0.01, 0.11) for deviations from R, H i — L 5100 and P H/3 ~ L H g relations, respectively. These 
correlations are shown in Figures 9 and 10. 
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For the Y — relation as shown by Figure 11, we obtained 

= 10.8i5 2 7 °, y sat = 3.0±g;| x 10 4 , and k 0 = 0.28 ± 0.03, (B4) 

from Figure 12, yielding the critical radius of 

r H/3 = 3.0+oi x 10 4 m° 7 - 29 ; or R cr it (H/3) = 17.1±|;| m\ 29 ltd. (B5) 

The results are very similar, and entirely consistent with the results of the average scheme presented 
in the main body of the paper. 
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8.— Mean and RMS spectra (observed flux vs. rest-frame wavelength) of the objects in the SEAMBH2013 
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Fig. 9. — Same as Figure 3 but for the direct scheme. The left panel and Figure 3a are identical, the right panel 
corresponds to Figure 3e. 
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Fig. 10. — Same as Figure 4 but for the direct scheme. The object with the largest error bar is NGC 7469. 
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Fig. 11. — Same as Figure 5 but for the direct scheme. 















































